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Summary. The study of gas and dust at high redshift gives an unbiased view 
of star formation in obscured objects as well as the chemical evolution history 
of galaxies. With today's millimeter and submillimeter instruments observers use 
gravitational lensing mostly as a tool to boost the sensitivity when observing distant 
objects. This is evident through the dominance of gravitationally lensed objects 
among those detected in CO rotational lines at 2: > 1. It is also evident in the 
use of lensing magnification by galaxy clusters in order to reach faint submm/mm 
continuum sources. There are, however, a few cases where millimeter lines have 
been directly involved in understanding lensing configurations. Future mm/submm 
instruments, such as the ALMA interferometer, will have both the sensitivity and 
the angular resolution to allow detailed observations of gravitational lenses. The 
almost constant sensitivity to dust emission over the redshift range z ~ 1 — 10 means 
that the likelihood for strong lensing of dust continuum sources is much higher than 
for optically selected sources. A large number of new strong lenses are therefore 
likely to be discovered with ALMA, allowing a direct assessment of cosmological 
parameters through lens statistics. Combined with an angular resolution < O'.'l, 
ALMA will also be efficient for probing the gravitational potential of galaxy clusters, 
where we will be able to study both the sources and the lenses themselves, free of 
obscuration and extinction corrections, derive rotation curves for the lenses, their 
orientation and, thus, greatly constrain lens models. 



1 Introduction 

Rapid progress in the development of millimeter astronomical facilities, such 
as the increase of antennae sizes and/or of the number of array elements, or 
as the continuing improvement in the sensitivity of detectors, have now made 
it possible to explore the high redshift universe in this window and therefore 
to exploit the potentialities of gravitational lensing effects. 

Why is it so important to explore this wavelength domain? In one short 
statement: the presence of cold dust and of molecular material can be traced 
in this window and both components witness the formation of heavy elements. 
If they are detected in galaxies at high redshifts, they allow us to probe star 
formation in the early universe. They reveal as well processes related to the 
startup of active galactic nuclei (AGN) and signal the presence of massive 
black-holes. Large amounts of molecular gas are encountered in the close 
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environment of the central engine in AGN. This material is often regarded as 
the fuel which allows to activate the AGN. An evolutionary scenario would 
then connect IR-luminous galaxies rich in molecular material and with intense 
star formation to the formation/feeding of massive black holes. 

Several fundamental questions are therefore underlying the search for 
dust and molecular gas at high redshift: the redshift of galaxy formation? 
the chemical evolution of the universe with time? the evolution of the dust 
content in the universe at early ages? the epoch and the scenario of the for- 
mation of massive black holes? the startup and evolution of AGN activity? 
Do we have already some clues to answer these questions? The most powerful 
AGNs, quasars, are now detected up to redshift around 6.5, and galaxies up 
to redshift 6. So we know that in this redshift range the universe already 
hosted galaxies and massive black holes and that its metal content was sub- 
stantial since the spectra of high redshift AGNs are very similar to those of 
low redshift AGNs. Yet, only a few objects are known at these high redshifts 
and this may provide a biased view. It is therefore mandatory to enlarge the 
sample and of course, the goal is also to push the redshift limit. 

Pushing the redshift limit also means that we are investigating sources 
with lower and lower flux density. This can be achieved through technical im- 
provements, using larger collectors and better detectors. The ALMA project 
is showing the way. Another manner is to take advantage of the effects in- 
duced by gravitational lensing (for a review of its theoretical basis, see the 
comprehensive book by Schneider et al. [126|). Firstly, image magnification 
allows us to detect more distant sources of cold dust and molecular gas of a 
given intrinsic luminosity, or to detect at a given redshift sources of fainter 
intrinsic luminosity. The latter in particular is important for good determina- 
tions of luminosity functions. Secondly, differential magnification effects can 
be used as an elegant tool to probe the size of molecular and dusty structures 
in the lensed source, as long as the lensing system provides the appropriate 
geometry. In this case, it is imperative to have an excellent model of the 
lensing system, as any structural information about the source itself for ex- 
ample is recovered by tracing the image back through the lensing system. 
Both aspects will be discussed at length in this paper. In some cases molecu- 
lar absorption lines allow us to obtain information about the lensing galaxy 
itself. 

Apart from hydrogen and helium, carbon and oxygen are the heavy ele- 
ments with highest abundance in the universe. Therefore, the CO molecule 
is the most suitable candidate for the detection of molecular gas in emission 
at high redshift. The CO molecules can be detected directly through their 
thermal line emission in the source or as silhouetted absorbers along the line 
of sight to a background source. The latter may occur for example for the 
lensing galaxy. Several other molecules have been detected at high redshift, 
HCO+, HCN, H2CO..., while dust is detected essentially through its thermal 
emission. 
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We provide in Sect. || an overview of the CO line emission and of the high 
redshift CO sources detected so far. The role played by gravitational lensing 
in studying CO sources at high redshift is highlighted. In Sect. ^ we review 
molecular absorption and the importance of such measurements to investigate 
the properties of the lensing galaxies. Sect. ^ introduces the dust continuum 
emission and the use of differential magnification effects which can be made 
to probe the dust content of the lensed objects. Three cases particularly well 
studied are discussed in detail in Sect. |5[ In Sect. ^ existing lens models for 
PKS1830-211 are reviewed and a new one introduced. Finally, in Sect. |^ the 
future of this type of investigations is presented in the perspective of new 
instrumental developments in general and of ALMA in particular. 

2 MOLECULAR EMISSION 

The goal of this section is primarily to highlight the benefits of exploiting 
gravitational lensing effects in the millimeter range, that is in CO line emis- 
sion. Therefore, after some brief comments on the pioneering observations 
of low redshift sources, we shall concentrate on the results obtained on high 
redshift sources. This section deals essentially with the CO line emission, 
while the following sections will discuss molecular absorption lines and dust 
thermal emission. 

Detection and measurements of the ^^CO rotational transitions in Galac- 
tic and extragalactic sources have had a great impact on the development 
of astrochemistry. The J=l-0 CO transition is excited by collisions with H2 
molecules, even in clouds at low kinetic temperature. The low-J transitions 
are in general optically thick, the opacity being determined observationally by 
the relative intensity of the corresponding ^'^CO transition. On the contrary, 
the high-J transitions are optically thin. Therefore, it is quite interesting to 
perform multi-transition studies to ascertain in a more secure fashion the 
physical conditions, temperature and density, in the emitting molecular ma- 
terial. This is particularly true in the case of very dense molecular material 
exposed to intense radiation fields, like in the environment of an AGN or in 
powerful star forming regions. The conversion factor which is used to derive 
the total mass of molecular material from the observed Leo, is also highly 
dependent on the physical conditions in the molecular material. It has been 
determined to be 4.6 Mq (K km s^^ pc^)^^ in standard Milky Way clouds 
[ p.38| . Recently it has been shown that lower values of the conversion factor 
(by up to a factor 10) should be used in the case of dense and warm material 
as usually encountered in a molecular torus around an AGN |Q . 

2.1 Low- and intermediate redshift galaxies 

Following far- infrared (FIR) observations by IRAS, an important population 
of IR-luminous (dust-rich) galaxies was found. This was the starting point for 
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investigating as well their molecular content, using the millimeter facilities 
available in the late 80's. A number of IR-luminous galaxies and AGNs were 
detected, mostly in the CO(l-O) transition and the field develop quickly. 
Regarding low redshift sources, let us briefly mention the detection of AGNs 
such as Mrk 231 |]l2|, Mrk 1014 p24|, I Zw 1 p and of some low redshift 
quasars [125| ||l| or radio galaxies ]114[ | |101| |100| . At moderate redshift, 
CO was detected in the radio galaxy 3C48 (z — 0.369) |127|. This search is 
continuing through e.g. the Caltech CO high and low redshift radio galaxy 
survey |t49| . 

On the side of CO sources at high redshift (z>2), the first object detected 
was IRAS 10214+4724 H [|36). All along the 90's, millimeter dishes and 
interferometer arrays in service were pushed to their limits in searching for 
other candidates, selected for example upon the strength of their submillime- 
ter flux. A large amount of observing time was dedicated to such programs at 
the OVRO, BIMA, Nobeyama and IRAM facilities. At face value, the success 
rate in detecting high redshift CO sources has been modest. One reason for 
this is the uncertainty in the precise redshift of the emitting molecular gas 
under search, while the backends of the instruments are narrow in compar- 
ison to the redshift range to be explored. Another reason is of course the 
limited sensitivity of existing instruments. Only the most luminous and most 
gas-rich systems can be detected. The situation should improve with new 
facilities such as ALMA (See Sect. 7.1). Sources detected so far are detailed 
below, in order of increasing redshift. 



2.2 High redshift galaxies 

The first source discovered, IRAS 10214+4724, at z = 2.285, has been de- 



tected in CO (3-2), CO (4-3) and CO (6-5) p§ [136|. A report on the detection 



of CO(l-O) [145| remains to be confirmed. The source is a gravitationally 
lensed ultraluminous IR-galaxy |Q The magnification factor is found 
to be around 10 for the CO source which has an intrinsic radius of 400 pc. 
Conversely, the far-IR emission detected in this object is magnified 13 times 
and arises from a source with radius 250 pc, while the mid-IR is magnified 
50 times and arises from a source with radius 40 pc. After correcting for 
magnification, and using a conversion factor L^q to M(H2) of 4 Mq (K km 



s pc ) [fl9|, the molecular gas mass is found to be 2 x 10 M©, in 



agreement with the estimated dynamical mass 3 x 10^° M©. As noted above 
however, such a value for the conversion factor, obtained from CO (1-0) ob- 
servations of Galactic molecular clouds, might not be applicable in the case 
of warmer and denser molecular material. A value for the conversion factor 
which is 5 times lower than the standard has been found in a study of ex- 
treme starbursts in IR-luminous galaxies |4^. Hence, the mass value quoted 
above should regarded as an upper limit to the mass of the molecular gas in 
IRAS 10214+472. Still pending is the question of the CO line emission share 
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Table 1. Galaxies at z>l with molecular emission (January 2002) 



Name 



Md/Mc 



-Jtir/L© Grav. lens Ref. 



BR1202-0725 4.69 

BR0952-0115 4.43 

BRI1335-0414 4.41 

PSS2322+1944 4.12 



6 X 10"' 2 X 10* 



3 X 10"* 3 X 10' 
1 X 10^^ continuum det. 
3 X 10" 1 X 10^ 



1 X 10' 



1 X 10' 



? 

YES 
NO? 
? 



log, |05 




APM08279+5255 


3.91 


2 X 10^ 


1 X 10^ 


~ 8 X 10'^ 


YES 




4C60.07 


3.79 


8 X 10"' 


2 X 10* 


~ 2 X 10'^ 


NO 




6C1909+722 


3.53 


4 X 10^° 


2 X 10* 


~ 2 X lO" 


NO 




MG0751+2716 


3.20 


8 X 10'° 






YES 






SMM02399-0136 


2.81 


8 X 10'" 


continuum dot. 


- 1 X 10'^ 


YES 


H 


MG0414+0534 


2.64 


5 X 10'° 


continuum det. 




YES 





SMM14011+0252 


2.56 


5 X 10'° 


continuum det. ■- 


- 3 X 10'^ 


YES 





H1413+117 


2.56 


2 X 10^ 


1 X 10* 


- 2 X 10'^ 


YES 





53W002 


2.39 


1 X 10'° 


weak continuum 




NO 


liH 


F10214+4724 


2.28 


2 X 10'° 


9 X 10* 


- 7 X 10'^ 


YES 






HR 10 



1.44 



7 X 10' 



2 X 10** 



' 9 X 10" NO 



Masses and fluxes corrected for gravitational magnification (approx.) 
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between a hidden AGN and a starburst in the 400 pc region surrounding the 
AGN. We notice also that the large extension (3 to 12 kpc) in CO emission 
reported by |l28t has not been confirmed by the IRAM interferometer data 
[|0|. 

One interesting case is that of the radio galaxy 53W002 aX z — 2.394, 
located at the center of a group of '^20 Lyman-a emitters. A possible detec- 



tion of the CO(l-O) line at Nobeyama was reported in [163|, al thou gh not yet 
confirmed by others. The first detection in CO(3-2) by OVRO [129| suggested 
a large extension (30 kpc) and the existence of a velocity gradient. None of 
these features has been confirmed by an IRAM interferometer data set with 
higher signal to noise ratio . From an astrometric analysis it is found that 
the 8.4 GHz and CO source are coincident, at a location consistent with that 
of the optical/UV continuum source. The most likely origin of the molecular 
emission is therefore from the close environment of the AGN. One should 
notice that 53W002 is definitely not a gravitationally lensed source. Using a 
conversion factor L^q to M(H2) in the range 0.4 - 0.8 M© (K km s~^ pc^)~^, 
more appropriate for dense and warm molecular gas around an AGN JTsf , 
the resultant molecular gas mass is found to be in the range (0.6 — 1.0) x 10^° 

Mq. 

The Cloverleaf, H1413+117 is a well known gravitationally lensed Broad 
Absorption Line (BAL) quasar aX z — 2.558 Bfl]. Its CO(3-2) transition was 



first observed with the IRAM 30m dish ||l|| and then with BIMA fieg. Later, 
the CO(4-3), CO(5-4) and CO(7-6) transitions have been detected, together 
with HCN(4-3) and a fine-structure line of CI. From a detailed analysis of 
these transitions, the molecular gas was found to be warm and dense 
with a low conversion factor of 0.4 M© (K km s~^ pc^)~^. High resolution 



maps in the CO(7-6) transition were obtained with OVRO |164| and with 
the IRAM interferometer |^ ||7^. The IRAM map has the best resolution 
and signal to noise ratio. Comparing with the HST images and exploiting 
differential gravitational effects, the CO (7-6) map allowed a derivation of 
both the the size and the kinematics of the molecular /d usty torus around the 



quasar central engine [[77| (for further details see Sect. After correcting 
for the amplification factor (30 according to the model used for the lensing 
system), and using the mean conversion factor 0.6 M© (K km s^^ pc^)~^ 
(derived by Barvainis et al. [^), the derived mass of molecular gas M(H2) 
is 2 X 10^ M©, in agreement with the dynamical mass of 8 x 10* M©. 

The source SMM 14011+0252 was observed with OVRO in CO(3-2) at 
z = 2.565 [p3, following its discovery as a strong submillimeter source de- 



tected in the course a survey of rich lensing clusters |134]. Correcting for an 
amplification factor of 2.75 and using a conversion factor of 4 M© (K km s~^ 
pc^)~^, the mass of molecular gas turns out to be 5 x 10^° M©, while the 
dynamical mass is found to be larger than 1.5 x 10^° M©. The CO emission 
is extended on scales of ^^10 kpc and associated with likewise extended radio 
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continuum emission [ [73[ . Optical and near- infrared (NIR) imaging shows two 
objects (designated Jl and 32), separated by 2'.'1 Q. Both the molecular gas 
and the radio continuum, however, have their strongest emission ~1" north 
of the J1/J2 components and are extended between Jl and J2. This suggests 
that the optical/NIR emission comes from two 'windows' in the obscuring 
molecular gas and dust and that J1/J2 represent emission from a coherent 
large galaxy. The extended nature of the radio continuum, the lack of X-ray 
emission and the lack of optical broad emission lines (Wiklind et al. 2002 
in prep) suggest that only star formation powers the large FIR luminosity. 
Assuming a Salpeter initial mass function and correcting for the gravitational 
magnification, the FIR luminosity indicates a star formation rate exceeding 
10^ Mq yr-i. 

The gravitationally lensed quasar MG 0414-1-0534 was observed with the 
IRAM interferometer in the CO(3-2) transition at z = 2.639 ||ll|. The lensed 
nature of this system is known from a 5 GHz map |Q. It displays four 
quasar-spots separated at most by ^ 2". The beam of the IRAM data (2'.'0 
X 0'.'9) does not allow to separate the components in the CO (3-2) velocity- 
integrated map. However, by fitting the UV data directly, it has been possible 
to resolve the combined A components (A1-I-A2) from component B and to 
get separate CO (3-2) spectra. The relative strength A:B in the 5 GHz radio 
continuum is 5:1. The millimeter continuum rather shows a ratio 7:1 and 
differences are seen between the A and B CO(3-2) spectra, suggesting that 
differential magnification effects may be at work. The magnification factor 
is unknown: hence only an upper limit can be derived for the molecular gas 
mass. Assuming in addition a conservative conversion factor of 4 Mq (K km 
s~^ pc2)~\ the upper limit found for M(H2) is 2.2 x lO" Mq. This figure is 
below the upper limit derived for the dynamical mass, 9 x 10^^ M©. 

Another source detected first through submillimeter observations is SMM 
02399-0136 |0. It is known to be gravitationally amplified by a foreground 
cluster of galaxies, the amplification factor being 2.5. This source has been 
detected with OVRO in the CO(3-2) transition at z = 2.808 The mass of 
molecular gas deduced in this object, correcting for a 2.5 amplification factor 
and using the conversion factor applicable to Galactic clouds, 4 Mq (K km 
s""'^ pc^)""'^, is 8 X lO"'^" Mq. From the upper limit on the apparent size of 
the CO emitting source (5") and the width of the CO line (710 km s~^), an 
upper limit to the mass of molecular gas of 1.5 x lO^*' Mq can be derived |]53| . 
The SCUBA results indicate that SMM 02399-0136 is an IR-hyperluminous 
galaxy. On the other hand, optical data shows clearly than it hosts as well a 
dust-enshrouded AGN |^ . A precise share of CO emission between the two 
components remains to be investigated. 

In the course of a systematic CO emission survey of gravitationally lensed 
sources with the IRAM interferometer ||l^, the source MG 0751-1-2716 has 
been detected in the CO(4-3) transition a.t z — 3.200. This source was first 
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discovered to be a gravitationally lensed quasar, from VLA maps [Q. It 
shows four quasar-spots with maximum separation of '.'9. The lensing galaxy, 
which provides the image geometrical configuration, is part of a group of 



galaxies adding another shear to the lens-system |144|. The lensing system 
remains to be modeled in detail. Therefore, the amplification factor is not 
known. However, given the observed strength of the CO line emission it should 
be large. Assuming that the CO emission in this source is mostly from the 
close environment of the AGN, we consider a conversion factor in the range 
0.4 to 1 Mq (K km s~^ pc^)~^. The corresponding upper limit (no correction 
applied for the unknown amplification factor) for the mass of molecular gas 
is in the range 8 x lO^" to 2 x 10" Mq @. 

The distant powerful radio galaxy 6C 1909+722 has been detected in the 
CO(4-3) line at z = 3.53, with the IRAM interferometer, and in dust sub- 



millimeter emission using SCUBA 1 107 . It is unlikely to be a gravitationally 



lensed object. Hence, the derived mass of molecular material is quite large, 
even assuming a conservative value for the conversion factor (about one fifth 
the value derived from Galactic molecular clouds). It is found to be in the 
range (0.5 - 1.0) x 10" M, 



0- 



Another possibly unlensed powerful radio galaxy 4C 60.07, has been de- 
tected in the CO(4-3) line emission at 2; = 3.79 at IRAM, and in dust thermal 
emission at submillimeter wavelengths with SCUBA and at millimeter wave- 



lengths at IRAM |107|. Remarkably, the CO line emission extends over 30 
kpc and breaks into two components: one corresponding to the AGN (radio 
core) and a second one which is rather related to a major episode of star 
formation. This state of merging is speculatively interpreted as the formative 
stage of an elliptical host around the residing AGN. Again, the molecular 
mass is found to be quite large, around 10^^ M0. 

The gravitationally lensed BAL quasar, APM 08279+5255, has been de- 
tected in the CO(4-3) and CO(9-8) transitions at z = 3.911, with the IRAM 
interferometer [14 3| . The CO line ratio points towards warm and dense molec- 
ular gas. Thermal emission from the dust component is also measured. Both 
the molecular and dust luminosities appear to be very high. Gravitational 
amplification is therefore suspected and has subsequently been confirmed 
through the detection of three optical/NIR components j85| (see also 
Sect, and Fig. |l3|). The magnification factors for the molecular gas and 
dust where estimated After correcting for these factors, the dust mass 
is found to be in the range (1 — 7) x 10^ M©, and the molecular gas mass 
in the range (1 — 6) x 10^ M©. In this interpretation, the molecular/dusty 
component is in the form of a nuclear disk with radius 90-270 pc orbiting the 
central engine of the BAL quasar. This source looks therefore quite similar 
to the Cloverleaf. Recently, however, extended low-excitation CO emission 



(the J=l-0 and J=2-l transitions) have been detected using the VLA |106| 
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This extended emission is likely to be associated with a cooler molecular 
component than the CO (9-8) emission. 

Finally, let's discuss the four sources detected so far at z larger than 4: 

PSS 2322+1944 The radio quiet quasar PSS 2322+1944 has recently been 
detected in the CO (5-4) and CO (4-3) transitions at a redshift oi z ~ 4+2 
with the IRAM interferometer The velocity- integrated CO line fluxes are 
3.74 ± 0.56 and 4.24 ± 0.33 Jy km s"\ with a linewidth « 330km s^^. The 
1.35 mm (250/im restwavelength) dust continuum flux density is 7.5 mJy, 
in agreement with previous measurements at 1.25 mm at the IRAM 30m 



telescope [|104| , and corresponds to a dust mass of ~ 10^ Mq. With the present 
angular resolution of the observations, no evidence for extended emission has 
been found yet. The implied gas mass is estimated to be « 3 x 10"^^ Mq, using 
a conversion factor of 4.6 Mq K km s^^ pc^. The properties of PSS 2322+1944 
are described in detail in |Q . 

BRI 1335-0415 BRI 1335-0415 was detected in the CO(5-4) transition with 
the IRAM interferometer at 2 = 4.407 |6^. The source does not exhibit a 
noticeable extension neither in the 1.35mm continuum nor in the CO line 
emission. In addition, there is no obvious sign of gravitational lensing on 
the line-of-sight to this source. The authors have derived a very large mass 
of molecular gas, close to 10^^ Mq. Even with a conversion factor 3 times 
smaller, more appropriate for this type of object, the mass remains a few lO^'' 

Mq. 

BR 0952-0115 The gravitationally lensed radio quiet quasar BR 0952-0115 



has been detected in CO (5-4) with IRAM facilities at z = 4.43 |6l|. A tenta- 
tive estimate of the mass of molecular material M(H2) is 3 x IO^Mq. Note 
however that a more precise model of the lens has still to be worked out. 

BR 1202-01215 BR 1202-01215 is the most distant source detected in CO. 
It has been reported in the CO (5-4) transition observed with the Nobeyama 



array flO^], and in the CO (4-3), CO (5-4) and CO (7-6) lines observed with 



IRAM facilities |105]. The CO maps show two separate sources on the sky: 
one is coincident with the optical quasar (for this source CO (5-4) provides 
z — 4.695), while the other is located 4" to the North- West, where no optical 
counterpart is found (for this source CO (5-4) provides z ~ 4.692). At this 
redshift the 4" extension corresponds to a de-projected distance of 12-30 
kpc. It is uncertain whether this object is gravitationally lensed or not. Is 
the North- West source a second image of the quasar? There are hints that 
this might be the case as a strong gravitational shear has been measured 
in the field (Fort and D'Odorico, private communication). Else, each of two 
separate sources ought to have its own heating source, AGN or starburst. 
Assuming no gravitational boost and using a conversion factor of 4 Mq (K 
km s^^ pc^)^^, the mass of molecular gas is quoted to be 6 x 10^° Mq. 
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A summary of the sources properties is provided in Table |l]. This quick 
compilation of high redshift CO sources has prompted a number of key-issues: 

(i) A major difficulty in detecting high redshift CO sources is the lack of 
precision in our guess for the redshift of the molecular gas. The instantaneous 
frequency coverage of current backends requires that the redshift is known 
a-priori with a precision of a few percent. Why is this condition hard to 
fulfill? The molecular gas emission in distant objects can arise from the close 
environment of an AGN. Yet, published redshifts for distant AGN are mostly 
measured from emission lines of highly ionized species which can be strongly 
affected by winds. Indeed velocity offsets of up to 2500 km have been 
observed between the CO lines and the blue-shifted CIV line for example 
(e.g. [Q). Hopefully, this limitation will be overcome with the next generation 
of backends. 

(a) An uncertain part in the interpretation of the observed CO line intensi- 
ties lies with the physical state of the molecular gas and the conversion factor 
L^Q to M(H2) to be applied in the case of high redshift sources. When sev- 
eral CO transitions are observed (such as for IRAS10214+4724, H1413+117, 
APM08279H-5255 and BR 1202-0725), the physical conditions of the molec- 
ular gas can be analyzed, pointing towards warm (T~100 K), dense (a few 
10^ cm~^) and moderately optically thick material. Such conditions could 
very well characterize molecular gas in the proximity of an AGN. Conversely, 
the conditions of the molecular gas in an extended starburst may be more 
similar to those encountered in Galactic molecular clouds. The conversion 
factor depends on the physical conditions of the molecular gas. It ranges 
from a value of 4 (K km s^^ pc^)^^ in Galactic clouds [I19|, to 1 
(K km s""'^ pc^)~^ in IR-ultraluminous galaxies | |4l|] and possibly 0.4 M© (K 
km pc^)""'^ in the surroundings of an AGN |13| . Therefore, it would be 
important to have some clues about the share AGN/starburst in the heating 
mechanism for high redshift CO sources. In that respect, the CO line width 
and the compactness of the source may bring some pieces of information. 

(in) Finally, the outmost efforts should be made to find out whether a source 
is gravitationally lensed or not, before the claim for the presence of a huge 
amount of molecular gas (10^^ Mq) can be taken as a starting point for 
modeling. Weak shear from an intervening galaxy cluster, like in the cases of 
SMM 14011+0252 and SMM 02399-0136, induces a mild magnification factor 
in the range 2-3. Strong shear (possibly combined with weak shear), induces 
magnification factors of up to 30! This would decrease by one order of mag- 
nitude the mass of molecular gas derived. If, at the same time, the applicable 
conversion factor is on the low side of its possible values range, a reduction 
of the actual molecular gas mass by another order of magnitude would apply. 
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In conclusion, it is important to search for other high redshift CO sources 

and, at the same time, to investigate carefully the nature of their environment/line- 
of-sight and the physical conditions in their molecular gas. 



Another method to study molecular gas at high redshift is to observe molec- 
ular rotational lines in absorption rather than emission. Whereas emission is 
biased in favor of warm and dense molecular gas, tracing regions of active 
high mass star formation, molecular absorption lines trace excitationally cold 
gas. This is important since a large part of the molecular gas mass may reside 
in regions far away from massive star formation and therefore remain largely 
unobserved in emission. 

Molecular absorption occurs whenever the line of sight to a background 
quasar passes through a sufficiently dense molecular cloud. In contrast to 
optical absorption lines seen towards most high redshift QSOs, molecular 
absorption is invariably associated with galaxies, either in the host galaxy of 
the continuum source or along the line of sight. In nearby galaxies molcciilar 
gas is strongly concentrated to the central regions, making the likelihood 
for absorption largest whenever the line of sight passes close to the center 
of an intervening galaxy. This, of course, means that molecular absorption 
in intervening galaxies is likely to be associated with gravitational lensing, 
and vice versa. Indeed, the only known systems of intervening absorption 
(B0218+357 and PKS1830-211) are gravitationally lensed and the absorption 
probes molecular gas in the lensing galaxy. Molecular absorption lines can 
thus be used to study the neutral and dense interstellar medium in lenses. At 
the present the sample of lens galaxies probed by molecular absorption lines is 
limited, but with the advent of a new sensitive interferometer instrument like 
ALMA, the number of potential candidate systems will increase substantially 
and make it possible to probe the molecular interstellar medium in the lens 
galaxies in some detail. Moreover, the molecular absorption lines provide 
unique kinematical information which is valuable when constructing a model 
of the lensed system. 

3.1 Detectability 

As mentioned above, molecular absorption traces a different gas component 
then emission lines. For optically thin emission the integrated signal Ico is 



where A^tot is the total column density of a given molecular species, 
the upper energy level of a transition with AE = hv, Tbg is the local 
temperature of the Cosmic Microwave Background Radiation (CMBR) and 
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J{T) = {tw/k){e'"'/'''^ ~ When T Tbg all the molecules reside in 

the ground rotational state J — and the signal disappears. For molecular 
absorption the observable is the velocity integrated opacity Ir^ : 



where TVtot is again the total column density of a given molecular species, 
while El is now the lower energy level and is the permanent dipole moment 
of the molecule. For the ground transitioij^, Ei — 0, Ir^ oc Ntot /i§ (1 — 
^-hu/kT^ Ri (hh'/k) Ntot fJ-oT~^- In contrast to emission lines, the observed 
integrated opacity increases as the temperature T decreases. 

Molecules are generally excited through collisions with molecular hydro- 
gen H2. The excitation temperature, Tx, therefore depends strongly on the 
H2 density. The collisional excitation is balanced by radiative decay and a 
steady-state situation with Tx = Tk requires a certain critical H2 density. For 
CO, which has a small permanent dipole moment /zq, the critical density is 
rather low, 4 x 10* cm~'^, while molecules with higher dipole moments require 
higher densities. For instance, HCO"^ which has a dipole moment more than 
30 times larger than that of CO, the critical density is 2 x 10^ cm""^. 

The strong dependence of the opacity on the permanent dipole moment 
means that absorption preferentially probes low excitation gas, i.e. a cold 
and/or diffuse molecular gas component. If multiple gas components are 
present in the line of sight, with equal column densities but characterized 
by different excitation temperatures, absorption will be most sensitive to the 
gas component with the lowest temperature. The dependence of the opac- 
ity on the permanent dipole moments also means that molecules much less 
abundant than CO can be as easily detectable. For instance, HCO"*" has an 
abundance which is of the order 5 x 10"'' that of CO, yet it is as easy, or 
easier, to detect in absorption as CO. This is illustrated in Fig. |l|, where the 
observed opacity of the CO(l-O) and HC0+(2-l) transitions a.t z ~ 0.25 are 
compared. In this particular case, the HCO^ line has a higher opacity than 
the CO line. 

3.2 Observables 

Analysis of the molecular absorption lines gives important information about 
both the physical and chemical properties of the interstellar medium. This 
can have implications for identifying the type of galaxy causing the absorption 
and, in some cases, help to identify the morphological type of lenses. In this 
section a short description of the analysis that can be done is presented. A 
more detailed description can be found in the references given in the text. 

^ This expression is strictly speaking only true for linear molecules. 
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Optical depth. The observed continuum temperature, Tc, away from an 
absorption line can be expressed as Tc = /sJ(T'b), where fs is the beam 
fiUing factor of the region emitting continuum radiation, Tb is the brightness 
temperature of the background source and J{T) = {hu/k)/[l — exp{—hiy/kT)] 
(e.g. [153|). The spatial extent of the region emitting continuum radiation 
at miUimeter wavelengths is unknown but is certain to be smaUer than at 
longer wavelengths. The BL Lac 3C446 has been observed with mm-VLBI 
and has a size < 30 /iarcseconds . Since the angular size of a single dish 
telescope beam at millimeter wavelengths is typically 10" — 25", the brightness 
temperature of the background source, Tb, is at least 10^ x Tc. This means 
that the local excitation temperature of the molecular gas is of no significance 
when deriving the opacity. The excitation does enter, however, when deriving 
column densities. 



Excitation temperature and column density. The excitation tempera- 
ture, Tx, relates the relative population of two energy levels of a molecule as: 
"2 _ S2 g-hv2i/kT„ where % is the statistical weight for level i and hvoi is 
the energy difference between two rotational levels. In order to derive T^ we 
must link the fractional population in level i to the total abundance. This is 
done by invoking the weak LTE-approximationQ. We can then use the parti- 
tion function Q{Tx) = X) jLo 9 J e"^-'/'^"^'' to express the total column density, 
iVtot, as 



9jAj,j+ 
Q{T^)e^-'/^^^ 



(1) 



where / T,ydV is the observed optical depth integrated over the line for a 
given transition, gj = 2 J + 3 for a transition J — > J +1, and Ej is the energy 
of the rotational level J. By taking the ratio of two observed transitions from 
the same molecule, the excitation temperature can be derived. The strong 
frequency dependence of the column density in Eq. is only apparent since 
the Einstein coefficient, Ajj+i, is proportional to f . 



3.3 Known Molecular Absorption Line Systems 

There are four known molecular absorption line systems at high redshift: 
z=0.25-0.89. These are listed in Table || together with data for the low redshift 
absorption system seen toward the radio core of Centaurus A. For the high 

^ In the weak LTE-approximation Tx ~ Trot, but the rotational temperature Trot 
is not necessarily equal to the kinetic temperature and can also be different for 
different molecular species. 
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Table 2. Properties of molecular absorption line systems. 



Source zi°' z<''' A^co Na^ N^^^ A;,'"*' Nui/Nh^ 

cm~^ cm"'^ cm~^ 



Con A 


0.00184 


0.0018 


1.0 


X 10^** 


2.0 X 


10"° 


1 X 


10"" 


50 


0.5 


PKS1413+357 


0.24671 


0.247 


2.3 


X lO^** 


4.6 X 


10^0 


1.3 > 


< 10"! 


2.0 


2.8 


B3 1504+377A 


0.67335 


0.673 


6.0 


X 10^** 


1.2 X 


10^1 


2.4 > 


< io"i 


5.0 


2.0 


B3 1504+377B 


0.67150 


0.673 


2.6 


X lO^** 


5.2 X 


10"° 


< 7 X 10"° 


<2 


<1.4 


B 0218+357 


0.68466 


0.94 


2.0 


X lO^^ 


4.0 X 




4.0 > 


< 10"° 


850 


1 X 10" 


PKS1830-211A 


0.88582 


2.507 


2.0 


X 10^** 


4.0 X 


10"" 


5.0 > 


< 10"° 


100 


1 X 10" 


PKS1830-211B 


0.88489 


2.507 


1.0 X 


IQie (■=) 


2.0 X 


10"° 


1.0 > 


< io"i 


1.8 


5.0 


PKS1830-211C 


0.19267 


2.507 


< 6 X lO^'^ 


< 1 X 


10"° 


2.5 X 10"° 


<0.2 


>2.5 



(a) Redshift of absorption line. 
(6) Redshift of background source. 

(c) 21cm HI data taken from |||] ||] [||| Q. A spin-temperature of 100 K and a 
area covering factor of 1 was assumed. 

(d) Extinction corrected for redshift using a Galactic extinction law. 

(e) Estimated from the HCO"*" column density of 1.3 x 10^'' cm"^. 



redshift systems, a total of 18 different molecules have been detected, in 32 dif- 
ferent transitions. This includes several isotopic species: C^^O, H^'^CO, 
H^'^CN and HC^^O''". As can be seen from Table |^, the inferred H2 column 
densities varies by ^ 10"^. The isotopic species are only detectable towards 
the systems with the highest column densities: B0218+357 and PKS1830-211, 
which are also the systems where the absorption originates in lensing galax- 
ies. The large dispersion in column densities is reflected in the large spread 
in optical extinction, ^v, as well as the atomic to molecular ratio. Systems 
with high extinction have 10-100 times higher molecular gas fraction than 
those of low extinction. 



Absorption in the host galaxy. Two of the four known molecular ab- 
sorption line systems arc situated within the host galaxy to the 'background' 
continuum source: PKS1413-fl35 Q and B3 1504-f 377 ||l5|. The latter 
exhibits two absorption line systems with similar rcdshifts, z=0. 67150 and 
0.67335. The separation in restframe velocity is 330 kms^^. This is the type 
of signature one would expect from absorption occurring in a galaxy acting 
as a gravitational lens, where the line of sight to the images penetrate the 
lensing galaxy on opposite sides of the galactic center. However, in this case, 
as well as for PKS1413+135, high angular resolution VLBI images show no 
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CO(l-O) 




Velocity (km/s' 



Velocity (km/s) 



Fig. 1. Plots of the observed opacity for the CO(l-O) and HCO'''(2-l) transitions 
seen at 2 = 0.25 towards PKS1413+135. The cut-off at opacities > 2 are due to 
saturation of the signals. The opacity of the HCO^(2-l) line is larger than that of 
the CO(l-O) line despite of an abundance which is 10"'^ — 10"^ that of CO. 



image multiplicity, despite impact parameters less than O'.'l (e.g. [113| |162|) 



The continuum source must therefore be situated within or very near the 
obscuring galaxy. 




Fig. 2. a) A 15 GHz radio image of the gravitational lens B0218-I-357 obtained 
with the VLA (courtesy A. Patnaik). b) A 15 GHz radio image of the gravitational 
lens PKS1830-211. 
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Fig. 3. Channel maps of HCO'''(2-l) and CO(3-2) absorption towards the 
B0218+357 obtained with the IRAM Plateau de Bure interferometer. The angular 
resolution does not resolve the two lensed images of the background QSO. The 
continuum weakens in the channel maps which corresponds to the absorption line, 
but never disappears completely. Since the absorption lines are strongly saturated, 
this shows that only part of the continuum is obscured by optically thick molecular 
gas. From Combes & Wiklind (unpublished). 



Absorption in gravitational lenses. The two absorption line systems 
with the highest column densities occur in galaxies which are truly interven- 
ing and each acts as a gravitational lens to the background source: B0218+357 
and PKS1830-211. In these two systems several isotopic species are detected 
as well as the main isotopic molecules, showing that the main lines are sat- 
urated and optically thick [Q |151| |l53|| . Nevertheless, the absorption 
lines do not reach the zero level. This can be explained by the continuum 
source being only partially covered by obscuring molecular gas, but that 
the obscured regions are covered by optically thick gas. The lensed images of 
B0218-I-357 and PKS1830-211 consist of two main components. By comparing 
the depths of the saturated lines with fluxes of the individual lensed compo- 
nents, as derived from long radio wavelength interferometer observations, the 
obscuration is found to cover only one of two main lensed components [150| 
[|l5l|. This has subsequently been verified through mm- wave interferometer 



data m |153| [140|. 



B0218+357 This is a flat-spectrum radio source lensed by an intervening 
galaxy. The lens nature was flrst identifled by Patnaik et al. |108|. The lens 
system consists of two components (A and B), separated by 335 milliarcsec- 
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onds (Fig ^). There is also a faint steep-spectrum radio ring, approximately 
centered on the B component. Absorption of neutral hydrogen has been de- 
tected at Zd = 0.685 showing that the lensing galaxy is gas rich. The 
redshift of the background radio source is tentatively determined from absorp- 
tion Hues of MgIlA2798 and H7, giving Zg w 0.94 |Q. Molecular absorption 
lines were detected in this system [150 further strengthening the suspicion 
that the lens is gas-rich and likely to be a spiral galaxy. The molecular ab- 
sorption lines do not reach zero level. Nevertheless, absorption of isotopic 
species show that the main isotopic transitions must be heavily saturated. In 
fact, both the ^^CO and C^^O transitions were found to be saturated as well, 
while the C^^O transition remained undetected |Q js^. This gives a lower 
limit to the CO column density which transforms to ~ 4 x 10^^ cm~^ 
and an Ay « 850 mag. 

That the molecular gas seen towards B0218+357 covers only one of the 
two lensed images of the background source can be seen in Fig. ||, where the 
continuum decreases at velocities corresponding to the absorption line but 
never completely disappears. Subsequent millimeter interferometry observa- 
tions have shown that the absorption occurs in front of the A-component, 
which is then expected to be completely invisible at optical wavelengths. 
Nevertheless, images obtained with the HST WFPC2 in broad V- and l- 
band, show both components (Fig. ^). While the intensity ratio A/B of the 
two lensed images is 3.6 at radio wavelengths [109|, A/B«0.12 at optical 
wavelengths. The V— I values show no significant difference in reddening for 
the A- and B-component. Hence, there is no indication of excess extinction 
in front of the A-component despite the large Ay inferred from the molecu- 
lar absorption. Since it is unlikely that the A/B intensity ratio is very much 
different at optical and radio wavelengths (differential magnification could 
introduce a small difference if the radio and optical emission comes from sep- 
arate regions) the A component appears sub-luminous in the optical. The 
other possibility is that the B component is over-luminous at optical wave- 
lengths by a factor 30 (or 1.4 magnitudes), possibly caused by microlensing. 
This latter explanation is, however, quite unlikely in view of the presence of 
large amounts of obscuring molecular gas in front of the A component. By 
compiling a sample of flat-spectrum radio sources from the literature, with 
properties similar to that of B0218-f357 (except the gravitational lensing 
aspect), correcting for different redshift and normalizing the observed lumi- 
nosities at 1/ = 10 GHz, it is possible to show that the optical luminosity 
of B0218-I-357 is abnormally weak |155| (Fig. ^). In this comparison the ob- 
served magnitude of the A component was used, multiplied by a factor 1.3 
in order to compensate for the B component using the magnification ratio of 
3.6. This clearly showed the A component to be sub-luminous, rather than 
the B component being over- luminous. The interpretation of this is that the 
A component is obscured by molecular gas, with an extinction that is very 
large. Some light 'leaks' out but through a line of sight which contains very 
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little obscuring gas, hence not showing much reddening in the V— I colors. 
Assuming that all the obscuration occurs in the A component, only ~3% of 
the photons expected from the A component reaches the observer. Since the 
extent of the optical emission region is very small, this suggests the presence 
of very small scale structure with a large density contrast in the molecular 
ISM of the lensing galaxy. 




Fig. 4. The gravitational lens B0218+357. Left : Optical image obtained with the 
HST in the I band (FW814), showing the A and B components (archival data). 
In contrast to the radio image (Fig. bp,), the A component is weaker than the B 
component. This is caused by obscuration of molecular gas, which gives rise to the 
observed molecular absorption lines at Zd = 0.688. Right : Normalized SED for 
flat-spectrum radio QSOs. The fluxes for the A component in B0218-I-357 as ob- 
served with the HST are marked by circles. Their positions suggest that only 3% of 
the flux supposed to come from this component reaches the observer. (From |155|). 



PKS1830-211 This is a radio source consisting of a flat-spectrum radio core 
and a steep-spectrum jet. It is gravitationally lensed by a galaxy at = 0.886 
[ |l5l| into two images of the core-jet morphology (Fig. ^d). The two cores are 
separated by 0'.'97 and the images of the jet form an elliptical ring. PKS1830- 
211 is situated close to the Galactic center and suffers considerable local 
extinction. Its lens nature was first suspected through radio interferometry 
[ |120| , but as neither redshift was known nor optical identification achieved 
(cf. its status as a gravitational lens remained unconfirmed. 

The lensing galaxy was found through the detection of several molecular 
absorption lines at Zd = 0.886 [151|. At millimeter wavelengths the flux from 
the steep-spectrum jets is completely negligible and it is only the cores that 
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contributes to the continuum. It was soon found that the molecular absorp- 
tion was seen only towards one of the cores, the SW image. However, weak 
molecular absorption was subsequently found also towards the NE image. 
This fortunate situation gives two sight lines through the lens and gives ve- 
locity information which can be used in the lens modeling (see Sect. ^). A 
second absorption hue system has been found towards PKS1830-211, seen 
as 21cm HI absorption at z — 0.19 |9^, making this a possible compound 
lens system. This intervening system complicates the lens models of this sys- 
tem. A potential candidate for the z = 0.19 absorption has been found in 
HST NICMOS images |7|. It is situated -4" SW of PKS1830-211 and is 
designated as G2. The molecular absorption lines towards PKS1830-211 and 
their use for deriving the differential time delay between the two cores will 
be described in more detail in Sect 5^ and Sect. ^. 



4 DUST CONTINUUM EMISSION 

The spectral shape of the far-infrared background suggests that approxi- 
mately half of the energy ever emitted by stars and AGNs has been absorbed 



by dust grains and then re-radiated at longer wavelengths |118 
The dust is heated to temperatures of 20-50 K and radiates as a modified 
black-body at far- infrared wavelengths. At the Ray leigh- Jeans part of the 
dust SED the observed continuum flux increases with redshift. This is known 
as a 'negative K-correction' and is effective until the peak of the dust SED 
is shifted beyond the observed wavelength range, which occurs at z > 10. 
Dust continuum emission from high redshift objects is therfore observable 
at millimeter and submillimeter wavelengths and is an important sources of 
information about galaxy formation and evolution in general and for gravi- 
tational lenses in particular. 

4.1 Dust emission 

Dust grains come in two basic varieties, carbon based and silicon based. Their 
size distribution ranges from tens of microns down to tens of Angstroms. The 
latter are known as PAH's (Polycyclic Aromatic Hydrocarbonates) . Except 
for the smallest grains, the dust is in approximate thermodynamical equi- 
librium with the ambient interstellar radiation field. The dust grains absorb 
the photon energy mainly in the UV and re-radiate this energy at infrared 
and far- infrared (FIR) wavelengths. The equivalent temperature of the dust 
grains amount to 15-100 K and they emit as an approximate blackbody. 

The spectral energy distribution (SED) of dust emission is usually rep- 
resented by a modified blackbody curve, F^, cx v'^B^{Td) (cf. ||142| JiSTC), 



where is the blackbody emission, Td the dust temperature and is the 
frequency dependence of the grain emissivity, which is in the range [3 — 1 — 2. 
Such representations have successfully been used for cold dust components 
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where a large part of the SED is optically thin. When t « 1 or larger, the 
observed dust emission needs to be described by the expression: 

F,^f2,B,iTd){l-e-^") , (2) 

where f2s is the solid angle of the source emissivity distribution, t^, is the 
opacity of the dust. Setting t^, = {v/vqY gives F^, oc i'^B^{Td) for ^ 1 
and Fi, oc B,^{T(i) for ^ 1. The critical frequency I'o is the frequency where 

= 1. 



The infrared luminosity. The total infrared luminosity is derived by 
integrating Eq. |2| over all frequencies. Here the flux density F^, corresponds 
to the energy emitted by dust only. The infrared luminosity for an object at 
a redshift z is given by 



LiR = 4^ (1 + zf Dl 



F^dv 



(3) 



where Da is the angular size distance^. The solid angle Q appearing in Eq. g 
is a parameter derived in the fitting procedure. In the event of a single dust 
component, Q can be estimated from the measured flux F^^ at a given rest- 
frame frequency Vr 



n 



6.782 X 10" 



GHz 



-3 



Jy 



1 



(4) 



Using some typical values (Td = 30 K, /3 = 1.5, = & THz (50^m), v,- = 
Vohs{^ + z) = 1400 GHz (j'obs = 350 GHz at z = 3) and, flnally, an observed 
flux of ImJy) we get J7 « 2 x lO"^**. For a spherical source with a radius 
r « Dx^JQIt:, this corresponds to a dust continuum emission region with 
an extent of only ~ 160 pc. 

Although this is a very rough estimate of the size of the emitting region, 
it shows, since typical observed values were used, that FIR dust emission 
from distant objects tend to come from very small regions. This will be of 
importance when considering the effects of gravitational lensing. 

^ Expressing Eq. ^ in a form directly accessible for integration, we get 



^ = 8.53 X 10" (1 



Da 
Mpc 



-(axf 



dx 



The integral can be integrated numerically with appropriate values of the param- 
eter a = kTd/huQ. 
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Fig. 5. The spectral energy distribution of two starburst galaxies, a) M82 and b) 
Arp220. Despite a difference in far-infrared luminosity of almost 2 orders of magni- 
tude, their spectral energy distribution are nearly identical. Notice also the presence 
of cold dust in M82, visible as an excess flux at millimeter and submillimeter wave- 
lengths. The SEDs have been fitted by a modified blackbody curve, which becomes 
optically thick at SOpm and which has /3 = 2.0 for M82 and f3 = 1.3 for Arp220 
and using a single temperature component of Td = 45 K for both galaxies. 



The dust mass. An estimate of the dust mass from the infrared flux 
requires either optically thin emission combined with a knowledge of the 
grain properties, or optically thick emission and a knowledge of the geometry 
of the emission region (cf. |6^). 

The grain properties are characterized through the macroscopic mass ab- 
sorption coefficient, k^. Several attempts to estimate the absolute value of 

as well as its frequency dependence have given different values (cf. p6|). 
Combining the same frequency dependence as in with a value given by 
[p4|, the mass absorption coefficient can be described as 



« 0.15 ( — — V'^ m^ kg~i 

V 375GHz / ^ 



(5) 



where corresponds to the restframe frequency. This expression corresponds 
to a grain composition similar to that found in the Milky Way. At frequencies 
where the emission is optically thin, the dust mass can now be determined 
fromR, 



Eq. ^ can also be expressed as: 



4.08 X 10* X 



Jy 



Da 



Mpc 



Dlil + zf. 



(6) 
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4.2 Detectability of dust emission 

A typical far-infrared spectral energy distribution (SED) of a starburst galaxy 
(M82) is shown in Fig. The SED of a more powerful starburst (Arp220) 
is shown in Fig. ||b. Perhaps the most striking aspect of these SEDs is their 
similarity, despite that they represent galaxies with widely different bolomet- 
ric luminosities. In both cases most of the bolometric luminosities comes out 
in the far-infrared: M82 has a far-infrared luminosity of 3 x 10^° Lq, while 
Arp220 is a so called Ultra-Luminous Infrared Galaxy (ULIRG) with a far- 
infrared luminosity of 1 x 10^^ Lq. The SEDs shown in Fig. |^ have been fitted 
by a modified blackbody curve, which becomes optically thick at SO/im and 
which has /3 = 2.0 for M82 and (3 = 1.3 for Arp220 (cf. Eq. |). The modified 
blackbody curves has been fitted using a single temperature component of 
T(j = 45 K for both galaxies. Notice, however, the presence of a colder dust 
component in the SED of M82, which is visible as an excess flux at millimeter 



and submillimeter wavelengths |143| 



The observed dust continuum emission originates from dust grains in dif- 
ferent environments and which are heated by different sources. Nevertheless, 
a remarkably large number of dust SEDs, like the ones shown in Fig. |^, can 
be well fitted by only one, or in some cases two dust components (cf. the cold 
dust component in M82). 

For a single dust temperature component, the flux ratio between the sub- 
millimeter (850/im) and the far-infrared (lOO/im) is strongly dependent on 
the dust temperature. For Td = 45 K (as in the case of M82 and Arp220), 
/ssoMm/Zioopm « 3 X 10"^, whilc for Td = 20 K, /sso^un/ZiooMm ~ 0.06, or 20 
times larger. Nevertheless, as long as the dust temperature is not extremely 
low, it is much harder to observe the long wavelength tail of the dust SED 
than the peak at ~100/im (except that in the latter case one needs to observe 
from a satellite due to our absorbing atmosphere). 

At millimeter and submillimeter wavelengths the SED can, to a first ap- 
proximation, be characterized by f^, oc v'^ , where 7 = 3 — 4. Hence, the 
observed fiux increases as an object is shifted to higher redshift. This effect 
is large enough to completely counteract the effect of distance dimming. An 
example of this is shown in Fig. ^, where the observed flux at 850/xm has 
been calculated for a FIR luminous, 5 x 10^^ L©, galaxy, for two different 
dust temperatures, Td = 30 K and Td = 60 K, and for two different cosmolo- 
gies. The largest uncertainty in the predicted ffux as a function of redshift 
comes from the assumed dust temperature, rather than the assumed cosmol- 
ogy. However, regardless of dust temperature and cosmology, the effect of the 
'negative K-correction' of the dust SED is to make the observed flux more or 
less constant between redshifts of z = 1 and all the way to z w 10, where the 
Wiener part of the modified blackbody curve is shifted into the submillimeter 
window and the flux drops dramatically. 

This constant flux over almost a decade of redshift range makes the mil- 
limeter and submillimeter window extremely valuable for studies of the for- 
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mation and evolution of the galaxy population at high redshift in general and 
for gravitational lensing in particular. For a constant co-moving volume den- 
sity, the submm is strongly biased towards detection of the highest redshift 
objects. The prerequisite is, of course, that galaxies containing dust exist at 
these large distances and that the low flux levels expected can be reached 
by our instruments. Both of these criteria are actually fulfilled; powerful new 
bolometer arrays working at millimeter (MAMBO, and recently SIMBA) and 
in the submillimeter (SCUBA) have shown that low flux levels can be ob- 
served and that obj ects containing large amounts of dust do exist at early 
epochs (cf. (67| |l3| ||). 




Redshift 

Fig. 6. The observed flux density at 850/im of a source with Lpir = 5 x 10^^ L© 
as a function of redshift. The top set of curves correspond to a dust temperature 
of 30 K and the bottom curves to a dust temperature of 60 K. The full drawn 
lines correspond to a flat matter dominated universe (i7m = 1.0, = 0), and the 
dashed curve to a flat A dominated universe = 0.3, Oa ~ 0.7). 



4.3 Submillimeter source counts 

One of the first studies using long wavelength radio continuum emission was 
to simply count the cumulative number of detected sources as a function of 
flux level. These observations mainly probed high luminosity radio galaxies 
and showed a significant departure from an Euclidean non-evolving popula- 



tion. This was the first evidence of cosmic evolution |122| | |74| . 

The negative K-correction in the mm-to-far-infrared wavelength regime 
for dust emission has enabled present day submm/mm telescopes, equipped 
with state-of-the-art bolometer arrays, to get a first estimate of the source 
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counts of FIR luminous sources at high redshift. There are two bolometer 
arrays which have produced interesting results so far; SCUBA on the JCMT 
in Hawaii and MAMBO on the IRAM 30m telescope in Spain. Additional 
arrays are under commissioning and will likely contribute to this area shortly: 
SIMBA on the 15m SEST on La Silla, and BOLOCAM on the 10m CSO on 
Mauna Kea. 

The SCUBA bolometer array at the JCMT was put to an ingenious use 
when it looked at blank areas of the s ky ch osen to be towards rich galaxy 
clusters at intermediate redshift [132| [134| |18 (l9). The gravitational 



magnification by the cluster enabled very low flux levels to be reached and 
several detections were reported. This method has been used by others as 
well and an example of an image of the rich cluster Abell 2125 at 1250/Lfm is 
shown in Fig. Q |Q. More than a dozen sources are detected above the noise 
but none is associated with the cluster itself. Instead they are all background 
sources gravitationally magnified by the cluster potential. 

The cumulative source count of a population of galaxies is simply the 
surface density of galaxies brighter than a given flux density limit. In a blank 
fleld observation it is in principle derived by dividing the number of sources 
with the surveyed area. The effects of clustering has to be considered if the 
observed area is small. In practice there are several statistical properties 
that have to be considered. Usually the threshold for source detection is not 
uniform across the mapped area. Since the sources are generally found close 
to the detector limit those which have fluxes boosted by spurious noise has a 
higher likelihood to be detected than those which experience a negative noise 
addition, which are likely to be lost from the statistics. This latter effect leads 
to an overestimate of the true source flux. The possibility of an additional 
bias through differential magniflcation will be discussed in Sect. ^7^ . 

The case of submm/mm detectedgalaxies behind foreground galaxy clus- 
ters is yet more complicated (cf. [Q). The gravitational lens distorts the 
background area and magnifles the source fluxes. The magnitude of these 
effects may vary across the observed field. A detailed mass model of the lens 
is needed in order to transform the observed number counts into real ones, as 
well as knowle dge about the redshift distribution of the sources. Small and 



collaborators ( |132| |133||134| ) initially observed 7 clusters, constructed or 
used existing mass models of the cluster potentials, and managed to obtain 
source counts at sub-mJy levels (cf. ) . Although the lensing effect of clus- 
ters allows observations of weaker fluxes, it introduces an extra uncertainty in 
the number counts. This is, however, not dominating the overall error budget 
[|8| . There is another beneficial effect with the lensing in that the extension 
of the background area alleviates the problem of source confusion. The angu- 
lar resolution of existing bolometer arrays is approximately 15" and source 
confusion is believed to be a problem at flux levels below 0.5 mJy. 

Other blank field surveys using SCUBA have pushed as deep as the cluster 
surveys, but without the extra magnification they probe somewhat higher flux 
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levels. Examples of such deep blank field surveys include the Hubble Deep 
Field North the fields used for the Canada- France Redshift Survey Q 
[ft5l, the Lockman hole and the Hawaii deep field region SSA13 |]. 

All these submm deep fields, including the cluster fields, are only a few 
square arcminutes. Using on-the-fly mapping techniques a few groups have 
recently started mapping larger areas but to a shallower depth (cf. |^ |131|)- 

Carilli et al. |Q combined the number counts from all the blank-field 
observations. The result is a cumulative source count stretching from ~15 
mJy to 0.25 mJy (Fig. ^). The source counts obtained using the lensing 
technique, after correcting for the lensing effects, are compatible with those 
obtained through pure blank-fields. The turnover at a flux level of ~10 mJy 
is probably real and represents a maximum luminosity of ~ 10^^ for an 
object at z « 3. The exact shape of the number counts is still uncertain at 
both the low and high flux ends. Results from the MAMBO bolometer array, 
which operates at 1250/im, have been multiplied by a factor 2.25 in order to 
transform it into the expected flux at 850/im. This assumes that the objects 
have an SED of the same type as starburst galaxies (cf. Fig. ||). 

In order to transform the cumulative source count into a volume density 
it is necessary to know the redshift distribution of the sources. It is, however, 
possible to circumvent this by fitting a model of galaxy evolution to the 
observed source counts. This has been explored extensively by Blain et al. p9[ , 



(see also 31 |141]), and will not be discussed further here. 



4.4 Submm source identification and redshift distribution 

The sources detected in submm/mm surveys can in a majority of cases be 
identified with sub-mJy radio sources (cf. ]135| ]). This population of weak 
radio continuum sources is believed to be powered by star formation rather 



than AGN activity |161| Attempts to identify the submm/mm sources 
with optical and/or infrared counterparts have failed in all but a small number 
of cases (cf. |^ |Q The submm/mm detected population is not related 

to nearby nor intermediate redshift sources, but are believed to be at 2 > 1, 
but the lack of clear optical/IR identifications has made it difficult to assess 
its true redshift distribution. An alternative technique for determining the 
redshift has been introduced by Carilli & Yun which relates the radio 
continuum fiux at 1.4 GHz with the measured flux at 850/im (see also 0). 
As the radio flux declines with increasing redshift, the submm flux increases 
(cf. Fig. H). Although the method is model dependent (mainly depending on 
the dust temperature Tdust, the radio spectral index as wel l as the frequency 
dependence of the dust emissivity coefficient, cf. Sect. [l.l[ ), it gives a rough 
estimate of the redshift. Using this method it has been possible to show that 
the majority of the submm/mm detected sources lie at a redshift 1 < z < 4 



(cf. 1 135 1 M) 



26 Wiklind & Alloin 




Fig. 7. An image of the cluster Abell 2125 obtained with the MAMBO bolometer 
array at the TRAM 30m telescope (from js^]). The angular size is in arcseconds. 
The noise rms is at 0.5 mjy/beam. 



4.5 Differential magnification 

One well-known property of gravitational lensing is that it is achromatic, 
meaning that the deflection of photons by a gravitational potential is inde- 
pendent of wavelength. The achromaticity is applicable to observed gravi- 
tational lenses as long as the source size is small compared to the caustic 
structure of the lens, such as when the Broad Line Region (BLR) of a QSO is 
lensed by a galaxy sized lens. Chromatic effects can, however, become impor- 
tant if the source is substantially extended (relative to the caustic structure) 
and the spectral energy density of the source is position dependent. 

The submm/mm detected dusty sources discussed in Sect. || are charac- 
terized by extended emission, several orders of magnitude larger than the 
compact sources generally studied in gravitational lensing. This applies to 
dust emission regardless whether the dust is heated by star formation or 
by a central AGN. Measured on galactic scales, however, the dust is rela- 
tively centrally concentrated, with typical scales ranging from 10^ pc to a 



few kpc (cf. Sect. 4.1). A dust distribution heated by a central AGN will 
have a radial dust temperature distribution, even when radiation transfer 
effects and a disk- or torus-like geometry are considered. This is observed 



in nearby Seyfert galaxies |115|. A radial temperature profile is also found 
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Fig. 8. Source counts from several surveys using SCUBA at 850^m and MAMBO 
at 1250/im (from [^). The dashed curve is a powerlaw of index —1.8 while the solid 
curve is an integrated Schechter luminosity function with a powerlaw index —2 and 
an exponential cut-off at lOmJy. All fluxes refer to 850/i. The 1250/im data points 
have been multiplied by a factor 2.25 in order to transform them into expected 
fluxes at 850/^m. 



in the case of a pure starburst |l30| , but spatially more extended than in 
the AGN case. Gravitational lensing of an extended dust distribution with a 
non-homogeneous temperature, and thus emissivity distribution, means that 
the assumption of achromaticity is no longer valid and the source may be 
differentially magnified. 

If the characteristic length scale in the source plane is 770 the characteristic 
length scale in the lens plane is ^0 = (Dd/Dg) rjQ. Taking a dust distribution 
of 1 kpc (770), a source redshift = 3 and a lens redshift — 1, the char- 
acteristic length scale in the lens plane becomes approximately 0'.'2. This is 
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close to the typical image separation for strong Icnsing. Since the submni/mni 
detected galaxies are believed to show a significant change in the dust tem- 
perature over this scale, it is quite likely that they will exhibit chromatic 
effects. 

An analytical model of the effect of differential magnification of dusty 
sources was presented in where it was shown that the effect can be 
strong and that it was most likely to produce an increase in the mid-infrared 
flux relative to the long wavelength flux. This would make the sources appear 
warmer than what their intrinsic SED would imply. 

A more detailed analysis of the effect of differenti al m agnification and 
its probability for occurance was done by Pontoppidan |116| . Using elliptical 
potentials and a realistic parameterization of the dust and its spectral energy 
distribution it was showed that both positive and negative distortions of the 
SED can occur. Here positive means an increase in the mid-IR part and 
negative means an increase in the far-IR/submm part. Fig. ^ shows a plot of 
the magnification in a cut through one of these models (elliptical potential) 
which does not hit the inner tangential caustic. By placing the center of 
a dust emission region, with a radial temperature profile, well outside the 
radial caustic (i.e. > |0'.'5| from the center), parts of the outer region of the 
dust distribution will fall on the high magnification plateau inside the radial 
caustic and be multiply imaged while the center is singly imaged and only 
moderately magnified. This situation would cause an enhancement of the 
long wavelength part of the SED relative to the mid-IR part. The radius of 
a typical dust distribution is ~ 0'.'2 at 2: « 3. If the center of the source is 
placed closer to the radial caustic, both the center and the extended dust 
distribution will be magnified, but the warmer central dust will experience a 
larger average magnification and hence result in a flattening of the SED at 
mid-IR wavelengths. Again, it might be that the cool dust is multiply imaged 
while the center (pos sibly containing an AGN) is singly imaged. 



Pontoppidan | 116 | found that the cross section for an enhancement of the 
long wavelength part of the SED is larger than for an enhancement of the 
mid-IR. However, the latter situation results in a stronger magnification and 
effects the observed SED to a higher degree. The latter case also represents a 
situation where the system is more likely to be recognized as a gravitational 
lens. 



Effect on number counts of submm/mm detected galaxies. A dis 

tortion map of the effect of differential magnification is shown in Fig 



10 



(from [116|). The caustic structure of an elliptical potential representing a 
central mass surface density of 3 x 10^ Mokpc^^, where the distortion of 
the dust SED due to differential magnification has been color coded. Black 
represents a negative distortion (cooler SED) and white represents a positive 
distortion (warmer SED). The dust distribution of the source is assumed to 
have a radius of 100 pc (left image) and 500 pc (right image). Quite naturally, 
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Fig. 9. Magnification in a cut through the caustic structure of an elliptical lens 
configuration. The two peaks corresponds to the radial caustic. This particular cut 
does not pass through the tangential caustic, which would have produced a yet 
stronger magnification peak close to the center. The length scale is in arcseconds. 
Notice the very strong gradient in the magnification when going from the one- image 
region to the the extended three- image 'plateau'. The magnification changes by a 
factor ~ 10 over angular scales of ~ O'.'Ol, corresponding to scales of ~ 50 pc at 
Zs « 2. From K. Pontoppidan's Master Thesis, Copenhagen University [116|. 



the larger the region over which dust is distributed, the larger is the region 
where negative distortion can occur. By placing the center of the source in 
the black/white regions of the distortion map, the observed SED will appear 
cooler/warmer. 

The implications for the submm/mm detected objects at high redshift 
is that differential magnification could induce a bias in the number counts. 
Especially since in most of the surveys done so far the sources are found close 
to the detection limit. The effect could induce an overestimate of the number 
of sources but it could also influence the slope of the cumulative number 
counts. The latter is more likely but better statistics from surveys reaching 
low noise levels are needed, as well as a better understanding of the total 
cross section for positive/negative distortions of dusty submm/mm sources. 

An interesting consequence of the differential magnification of these sources 
is that some, perhaps several, of the detected submm/mm objects may be 
multiply imaged systems when viewed at high angular resolution in submm/mm 
wavelengths. The radio identifications that have been done typically reach a 
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Fig. 10. The caustic structure of an elliptical potential showing areas where the 
dust spectral energy distribution will be influenced by differential magnification. 
The dust emission region is modeled as a circular disk with a radial temperature 
profile. The central heating can be either an AGN or a dense starburst. If the center 
of the dust emission region is placed in the black area, the observed SED will appear 
cooler than the intrinsic one, while the opposite effect occurs if the center is located 
in the black areas. The effect of a small (left) and large (right) dust region is 



illustrated. From K. Pontoppidan's Master Thesis, Copenhagen University |116 



resolution of 1" which is not sufficient to see multiple images on the expected 
0'.'l-0'.'2 scale. High angular resolution deep imaging with future instruments 
such as ALMA will resolve this issue. 



5 CASE STUDIES 

In order to describe in more detail the characteristics of millimeter obser- 
vations and interpretations of gravitationally lensed sources, as well as to 
illustrate their use, three cases are presented below. First is the luminous 
Broad- Absorption-Line (BAL) quasar APM08279-h5255 at z = 3.9. The grav- 
itational lens hypothesis for this source was put forward based only on its 
apparent luminosity. The second case is a detailed study of the quadruply 
lensed Cloverleaf quasar, where the gravitational lensing of molecular gas 
has enabled a more detailed and constrained lens model. The last example is 
PKS1830-211, where the lensing galaxy was actually first detected through 
millimctric molecular absorption lines at Zd = 0.886. The molecular absorp- 
tion lines in this system has been used to constrain the lens model by giving 
the velocity dispersion and are used to derive the differential time delay be- 
tween the two main lensed components. 
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5.1 APM08279+5255: A case of differential magnification? 

This object was discovered serendipitously during a search for Galactic car- 
bon stars l?^]. It was found to be a BAL QSO at a redshift z — 3.911 
(see for the redshift determination). With an astounding R-band mag- 
nitude of 15.2 and detection in three of the four IRAS bands, its bolometric 
luminosity turns out to be 5 x 10"'^^ L0. This in itself led to the suspicion that 
it is a gravitationally lensed object Subsequent observations, both from the 
ground and from space [|5| [|9|, led to the detection of three components, 
with a maximum separation of 0'.'35 ± 0'.'02, and with a flux ratio of the two 
brightest components of 1.21 ± 0.25 (cf. Fig. ^3|). The optical spectra of the 
two main components are similar to each other [8^ . No lensing galaxy has 
been identified, although the weak third image could potentially be the lens 
(see below). Nevertheless, based on the small separation of the main com- 
ponents, their similar spectra and the enormous luminosity inferred for the 
system, the lensing nature of this system is not questioned. Even in the case of 
strong gravitational magnification, APM08279-f 5255 is an intrinsically very 
luminous system, with Lboi ^ 10^'^ Lq. 

The high apparent brightness of APM08279-f-5255 has allowed a very good 
S /N optical spectra of the intervening absorption line systems to be obtained 
with the HIRES spectrograph on Keck Q . Several potential lens candidates 
are found as Mgll absorption line systems, with the most conspicuous one 
at z = 1.181. Placing the third image at this redshift, however, requires the 
lens to be unusually compact and luminous. It would need to be almost 5 
magnitudes brighter than an L* galaxy with the relevant velocity dispersion 
of ~ 150 kms"^ ||§ @- The possibility that the lens harbors an AGN can 
be dismissed since no emission lines from z < 3.9 are detected in the spec- 
trum. Also, the continuum of an intervening QSO should have been detected 
in the saturated parts of the absorption lines seen towards the background 
source. No such emission is detected (cf. @). APM08279+5255 could thus 
represent a 'text book' example of a gravitational lens with an odd number 
of components. 

Apart from being luminous at optical and UV wavelengths, APM08279+5255 
also contains large amounts of dust and metal rich molecular gas (Fig. ITl]) . 
The SED of APM08279-f 5255 is actually dominated by a strong dust con- 



tinuum emission (Fig. 12), detected over a wide wavelength band: from the 
restframe submm to mid-infrared bands. This puts APM08279-t-5255 in the 
class of hyperluminous IR galaxies even when correcting for a strong gravi- 
tational magnification. 

The overall dust spectral energy distribution is characterized by a steeply 
rising long wavelength part, with a change of slope around Arost = 200/im, 
and a flat mid-IR part. The dust continuum spectra can be fitted by two dust 
components. One 'cool' characterized by a dust temperature of Td = 200 K, 
which is optically thin at A > 200/im (cf. Fig. [T^ ). The second component 
is hot, with Td « 910 K, close to the sublimation temperature of carbon 
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based dust grains. This second component is optically thick. The total dust 
mass, uncorrected for gravitational magnification, is 2 x fO^ M0, most of it 
contained in the cool dust component. 

The CO emission lines shown in Fig. |ri| includes the high excitation tran- 
sition J = 9 - 8. The CO J = 9 level is J( J + 1) x 2.77 = 249 K above the 
ground state. Normal type Galactic molecular clouds with typical H2 densi- 
ties of ~ 300 — 10'^ cm"-^ are not sufficient to coUisionally populate the CO 
J = 9 level. The mere detection of the CO (9-8) line therefore shows that 
the gas has to be unusually dense and warm. This immediately suggests that 
this gas component resides close to the QSO, possibly associated with the hot 
dust component. If both the CO J = 4 — 3 and J = 9 — 8 emission are asso- 
ciated with the same gas component, the total molecular gas mass, corrected 
for magnification, is quite modest: 3 x 10^ M0 Q. If the lower transition, 
on the other hand, emanates from a more extended and cooler region than 
the J = 9 — 8 transition, the total molecular gas mass can be one to two 
orders of magnitude larger. That this is likely to be the case was shown by 
the detection of CO J = 1 - and J = 2 - 1 emission from APM08279-H5255 
(Papadopolous et al. 2001)^. Using the same conversion factor between II2 
column density and velocity integrated CO intensity as is used for the Milky 
Way and nearby galaxies, the total molecular gas mass in APM08279-f 5255, 
uncorrected for gravitational magnification, is (0.6 — 3.2) x 10^-'^ M© ]106| |. The 
amount of gravitational magnification is in this case expected to be low due to 
the extended nature of the molecular gas, especially the gas seen in the lower 
transitions. Incidentally, three additional CO emitting sources are detected 
within 3" of the center of APM08279-h5255 [[06|. If these are not gravita- 
tionally lensed images, which they are not if the currently best lens models 
are used (cf. ||6^ |Q ) , these three additional sources are not magnified to any 
significant degree. The field around APM08279-I-5255 should then represent 
a remarkable over-density of gas rich galaxies at high redshift. These three 
additional sources are, however, not detected in continuum emission with the 
Plateau de Bure interferometer nor at optical or NIR wavelengths and their 
exact nature remains undetermined. 

APM08279-I-5255 has a SED which is essentially flat from a restframe 
wavelength of ~30/xm to optical wavelengths (cf. Figjl^). This is usually 
interpreted as being the effect of a face-on configuration of a dust-disk sur- 
rounding a central AGN. The low inclination of the disk enables the observer 
to get an un-obscured view of the hot dust close to the AGN as well as the cool 
dust further away. Comparison with dust models calculated by Granato et al. 
[p8| [p9[ shows that the mid-IR slope is too shallow even for the most extreme 



^ The offset between the CO emission presented in Papadopoulos et al. |106| and 
that of Downes et al. results from the use of slightly different coordinates for 
APM08279-I-5255. The coordinates given in the caption of Fig. corresponds to 
the best optical/IR coordinates determined from both ground and space based 
imaging. 
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face-on models |89|, i.e. the dust SED in APM08279+5255 appears to be too 
'warm' even if heated by a powerful AGN. Another possible explanation for 
the flat mid-IR SED is that APM08279-I-5255 experiences differential magni- 
fication of the dust emission region. This possibility was explored by Egami 



et al. |47[ by applying sources of various sizes to their lens model. In Fig. 12 



the SED of APM08279-h5255 is shown together with a starburst model |121[ . 
The starburst model has been arbitrarily fitted to the long wavelength part 
of the observed SED. At mid-IR wavelengths, the starburst model predicts a 
flux which is ~50 times lower than the observed fluxes in APM08279-I-5255. 
The data points marked by open circles are the equivalent observed fluxes 
diminished by a factor of 50 in order to fit the starburst model. Although 
APM08279-I-5255 undeniably contains a powerful AGN, which is likely to 
contribute a substantial part of the heating of the gas and dust, the influ- 
ence of star formation can not be ruled out. Can differential magnification 
account for at least part of the difference between a pure starburst SED and 
the observed one? 




Fig. 11. CO spectra and maps of APM08279-I-5255 observed with the IRAM 
Plateau de Bure interferometer (Downes et al. [^^). Upper left: 1.4mm dust con- 
tinuum. Lower left: CO(9-8) emission. Upper and lower right: CO(4-3) and 
CO(9-8) emission line profiles. The angular resolution of the maps is 3'.'2 x 2'.'3, far 
too coarse to resolve the individual components seen at optical/NIR wavelengths. 
The maps are centered on 08''31'"4P!70, -F52°45'17'.'35 (J2000). 
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Fig. 12. The spectral energy distribution of APM08279+5255 fitted by two isotlier- 
mal greybody models. The full-drawn line corresponds to the sum of a 'cool' com- 
ponent (Tduat = 200 K) and a hot component (Tdust = 910 K). The dust emission 
becomes optically thick at A < 200^m. The data points (filled circles) are from 
Irwin et al. [Q, Ledoux et al. |^], Lewis et al. jsoj, Downes et al. |^^, Egami et 
al. ^1 and the Faint Source IRAS catal og (N ED). The dotted line is a starburst 
model of Rowan- Robinson & Efstathiou [121|, which has been arbitrarily fitted to 
the long wavelength part of the SED. The data points marked by open circles are 
the equivalent observed data points but with their values reduced by a factor 50 in 
order to fit on the starburst model. 



Modeling the lens APM08279+5255. The lensing configuration of this 
system has been modeled by Egami et al. and Ibata et al. |69|. Using an 
isothermal elliptical potential with no external shear, two different types of 
lens models were applied: a three-image model and a two-image model. The 
former model is non-singular in order to produce the third image, while the 
latter assumes that the third image is the lensing galaxy and the potential 
is singular in order to suppress the formation of the third image. The two- 
image model produce a modest magnification of ^7 (cf. |^), while the three- 
image model produce a magnification of ^90 for both a point source and a 
more extended source distribution |^ Since the apparent bolometric 
luminosity exceeds 10^^ Lq, the three-image configuration is more appealing. 
However, the core radius is large, 0'.'21, almost as large as the Einstein radius, 
0'.'29. If the lens is at a redshift k, 1.2, the core radius corresponds to ~1.2 
kpc in the lens. This is much larger than most measured core radii. In a 
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Fig. 13. The 2-image and 3-image lens model solutions for APM0927+5255 and 
their effect on the spectral energy distribution. In the 2-image solution (top left), 
the intrinsic mid-IR part of the SED is slightly depressed, while in the 3-image 
solution (bottom left) the intrinsic mid-IR is strongly enhanced (representing a 



positive distortion, see Sect. i.5). The intrinsic SED is represented by dashed and 



dotted lines, for a dust distribution with a radius of 600 pc and 300 pc, respectively. 
The full drawn line is the observed SED. In the top right is APM08279-I-5255 seen 
with the HST WFPC2 camera and lower right is the dust model and its location 
used in the 3-im age solution. From K. Pontoppidan's Master Thesis, Copenhagen 
University 



survey of 42 giant elliptical galaxies it was found that the objects which can 
be resolved have a median core radius of 225h^^ pc [Q. In the case of a 
three-image model, the potential is almost circular with e — 0.012, while the 
two-image model gives e = 0.083 The difference in ellipticity corresponds 
to a difference in the size of the caustic structure, which in turn influences 
the effects of differential magnification. In the three-image model, the caustic 
structure is approximately 45 pc in extent in the source plane, while the 
two-image model has a caustic structure almost 5 times larger. 

The effects on the lensing behavior for a source with a finite extent was 
explored in [M. A source with an extent exceeding ~ 500 pc resulted in 



a filled disk. A more detailed model was done by Pontoppidan |116] where 
an assumed source temperature distribution was used in order to derive the 
resulting spectral energy distribution. Using the model parameters of pTf, 
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where the source is located between the radial and tangential caustic for the 
three-image model, the hot dust is expected to be moderately enhanced by 
the outer magnification plateau (cf. Fig. In the two- image scenario, the 
QSO is again located outside the tangential caustic. In this case, however, 
the radial caustic is lacking due to the singular potential. The latter scenario 
can produce a modest negative distortion of the SED (as seen in the top left 
panel of Fig. |l^). In the three-image scenario, however, the effect on the SED 
is more dramatic and represents a positive distortion, i.e. the mid-IR part 
of the SED is enhanced relative to the long wavelength part (bottom left 
panel of Fig. |l^). The magnitude of the distortion is quite large, its details 
depending on the extent of the dust region. For a dust distribution with a 
radius of 650 pc, the differential magnification can enhance the intrinsic flux 
at restframe mid-IR wavelength with a factor ~ 10. A smaller extent of the 
dust results in a smaller enhancement factor. The dust region (for a radius of 
300 pc) and the caustic structure are seen in the lower right panel of Fig. 

The three-image model of APM08279-I-5255 is more likely to be correct 
than the two-image model since it produces a magnification which corre- 
sponds to a source with a bolometric luminosity which is large, but not ex- 
treme. The three-image model also means that for realistic dust distribution, 
the restframe mid-IR is strongly enhanced relative to the longer wavelength 
part of the SED. The intrinsic SED of APM08279-f 5255 resembles that of 
less extreme dusty QSO spectra (cf. [|o|). In fact, the shape of the intrinsic 
SED of APM08279-I-5255 now resembles that of pure starburst models, ex- 
cept that the mid-IR is still enhanced by a factor 5 — 10, marking the influence 
of the AGN. This shows that the effects of differential magnification must be 
considered before applying radiation transfer models to gravitationally lensed 
dusty sources. 

In order to model the resolved and extended low-J CO emission, the 
effects of a highly elliptical lensing potential has been explored The 
lensing galaxy is here assumed to be an edge-on spiral. The result is a good 
fit with the extended low-excitation CO emission, while the point sources 
from the background QSO, although imaged into three components, have 
widely different magnification ratios compared to the observed values. This 
may not be of great importance if microlensing affects the optical photometric 
results (e.g. M). 



5.2 The Cloverleaf: Another case of differential magnification 

The Cloverleaf is the gravitationally lensed image of the BAL quasar H1413-I-117 
at z = 2.558, showing four quasar-images (hereafter called spots) with angu- 
lar separation from 0'.'77 to 1'.'36. Since its discovery |Q, the Cloverleaf has 
been imaged with ground based telescopes in numerous bands up to I and 
with HST/WFPC2 in the UV, optical and near-IR [|46) Q |7|. 
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The lensing system. After the early lens model of Kayser et al. these 
new data sets have been used to derive an improved model of the lensing 
system Q virhich now includes: 

1. A cluster of galaxies with derived photometric redshifts in the range 0.8 
to 1.0, which contributes to the magnification. 

2. A lensing galaxy close to the line of sight to the quasar, which determines 
the geometry of the image (four main spots) and carries the largest share 
of the magnification. The redshift of the lensing galaxy has been ten- 
tatively measured with VLT/ISAAC at a value of 0.9 (Faure et al. in 
prep). 

In the following, we use this new model for the lensing system, which is 
essentially constrained by the HST data. Further details can be found in 



The IRAM milUmeter data sets. After its discovery in the CO (3-2) 
line emission with the IRAM Pico Veleta dish |12l , the Cloverleaf has been 



observed in the millimeter range by various teams and instruments ([160 



with BIMA, ||16J with OVRO). Yet, the best data sets collected to date 
on this object are from the IRAM Pico Veleta dish and Plateau de Bure 
interferometer. 

A total of six millimeter transitions have been reported from observa- 
tions with the IRAM Pico Veleta dish: CO (3-2), CO (4-3), CO (5-4), CO (7-6), 
CI(3Pi-3po) and HCN(4-3) |l|. Detailed non-LTE modeling of the CO line 
strengths by these authors indicates that the molecular gas is warm (T larger 
than 100 K), dense ( n(Il2) density larger than 3 x 10'^ cm"^) and not very 
optically thick {tco < 3). These results suggest that the molecular material 
is close to a powerful heating source and might therefore be related to the 
environment of the central engine in the quasar. They also prompt us for 
not using the conventional conversion factor CO to II2 which is derived for 
molecular clouds in the disk of our Galaxy. 

Thanks to the strength of the CO(7-6) transition, a high resolution (0'.'5) 
map was obtained with the IRAM Plateau de Bure interferometer. A first 
CO (7-6) interferometric data set |^ has later been complemented with ob- 
servations at intermediate baselines |7^. The combined data has lead to the 
CLEANed map restored with an 0'.'5 circular beam, shown in Fig. |lj . In 
order to search for a velocity gradient, we have derived the spatially inte- 
grated line profile, following the procedure described in jSj. The CO (7-6) line 
profile (Fig. |l^) shows a marked asymmetry with a steep rise and excess of 
emission (with respect to a standard Gaussian) on its blue side and a slower 
decrease on its red side. Excluding the central velocity channel (so that the 
split in velocity is symmetric), we have built the blue (-225, -25 km/s) and 
the red {+25, +225 km/s) maps displayed in Fig. p^,d respectively. The 
difference between the red and blue CLEANed maps (Fig. p^) establishes 
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firmly the presence of a velocity gradient at the 8a level. Measurements of 
the characteristics of the spots from the CO (7-6) image have been performed 
(spot flux ratios, sizes and orientations) through a fitting procedure in the 
visibility domain, as explained in j^] . Final parameters are provided in Table 
1 of i fTZt where the spot sizes are intrinsic to the image, i.e. deconvolved by 
the interferometer beam: spots A, B and C definitely appear elongated. 




Fig. 14. Image of the Cloverleaf obtained with the IRAM Plateau de Bure interfer- 
ometer [Q. a) is the total CLEANed image, c) and d) are the blue and red part of 
the total emission profile, while b) shows the difference between the CLEANed blue 
and red images. A velocity gradient in the underlying source can be inferred from 
the residual. The data has been restored with a circular beam of angular resolution 
O'.'S. 



Comparing images in the UV and the millimeter range. Images 
in the UV/optical correspond in the quasar restframe to the emission from 
the accretion disk surrounding the quasar central engine. This latter source 
is expected to be point-like. The four spots on the HST images do indeed 
have a stellar- like appearance, being circular with a FWHM of about 0'.'068 
[ [77| . The absolute photometry and relative intensity ratios of the four spots 
(Table 4 in ) have been computed using the Sextractor software [Q . The 
large variation of the intensity ratios in U, compared to V, R and I bands can 
probably be explained by absorption along the line of sight by intervening 



Gravitational Lensing at Millimeter Wavelengths 



39 




Fig. 15. The ^^CO(7-6) spectrum from the z = 2.56 Cloverleaf quasar obtained 
with the IRAM Plateau de Bure interferometer . The thin line represents a best 
fit Gaussian profile. Notice the asymmetric line profile with excess emission at the 
blue part of the spectrum. 



galaxies. Alternatively, this effect can be ascribed to dust extinction at the 
redshift of the quasar, using an SMC-like dust extinction law [146|. The 
presence of such an absorbing medium in the close environment of the central 
engine could be put in relation with the BAL appearance of this quasar. 

The millimeter CO lines are expected to arise from an extended structure, 
the so-called dusty /molecular torus, with an intrinsic radius of a few 10 pc to 
a few 100 pc (according to models). In such a configuration, different parts of 
the extended torus will be positioned differently with respect to the caustic 
(the curve which represents in the source plane the signature of the lensing 
system). As the image properties and the amplification factor in particular, 
are ruled by the relative positioning of the source/caustic, the four spots on 
the CO(7-6) image, each corresponding to the extended torus, will be dis- 
torted with respect to the four spots on the HST image (corresponding each 
to a point- like source). This features what is called 'differential magnification 
effects' (see Sect. 4.5). From the blue bump appearing on the CO(7-6) line 
profile (Fig. |l^), we clearly see that the blue-shifted part of the CO line arises 
from a region of the molecular torus which is positioned closer to the caustic 
than the region emitting the red-shifted side of the CO line. In this way, we 
are able to recover detailed structural and kinematical information about the 
molecular torus in the quasar. 

In order to derive precisely the shear induced by the lensing system on 
the extended source in the quasar, it is imperative to register with a high 
accuracy the Cloverleaf image in a waveband corresponding to a point-like 
source in the quasar (accretion disk: UV restframe, that is an R band image 
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for example) and in a waveband corresponding to an extended source in the 
quasar (molecular torus: CO(7-6) line). The high precision required, better 
than 0'.'2, was achieved using a combination of the HST data and of CFHT 
data acquired over a larger field of view under extremely image quality . 




Fig. 16. The caustic structure and the CO source distribution in the Cloverleaf 
quasar |7^. The central ellipses represent the CO source distribution for the red 
and blue part of the emission profile, respectively. The dot in the center corresponds 
to the quasar UV point source. The scale on the axis is arcseconds. The result has 
been obtained by combining an HST optical image with the interferometric CO 
data. 



Derived properties of the molecular torus in the Cloverleaf BAL 
quasar at z — 2.558. We have used the model of the lensing system as 
constrained by the HST data and presented above. The total amplification 
factor for the CO emission is found to be 30. This amplification factor trans- 
late to a molecular gas mass M(H2) = 2 x 10^ Mq and an atomic hydrogen 
gas mass M(HI)= 2 x 10^ Mq jl3| 0. 

We have derived the properties of the molecular torus in the quasar us- 
ing the CO(7-6) maps: firstly using the total line flux and secondly, using 
separately each of the blue and red halves of the CO (7-6) line |7^. We find 
a typical size for the molecular torus of 150 pc (assuming Hg = 50 km s^^, 
fim = 1 and — 0). When we treat separately the maps corresponding 
to the blue-half line and red-half line, we find that the quasar point-like UV 
source is almost exactly centered between the region emitting the blue-half 
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of the CO line and the region emitting the red-half of the hne (Fig. 16). This 
is reminiscent of a disk- or ring-like structure orbiting the quasar at about 
75 pc and with a Keplerian velocity of 100 kms~^ (assuming a disk orienta- 
tion perpendicular to the plane of the sky). The resulting central dynamical 
mass would be about 8 x 10* M©. This value is in good agreement with the 
estimate of the molecular gas mass made above from the total CO line flux, 
provided uncertainties in the inclination of the molecular torus and in the 
conversion factor from I(CO) to N(H2). 

In conclusion, we can regard the case of the Cloverleaf as a first and 
enlightening example of what will become routine when ALMA becomes 
available. Indeed, exploiting differential magnification effects is an extremely 
promising technique. The effective angular resolution on the CO source in the 
quasar at z = 2.558, using this procedure, is ~0'.'03, or about 17 times smaller 
than the synthesized beam of the IRAM observations! And the amplification 
factor in this case is around 30! 



5.3 PKS1830-211: Time delay and the Hubble constant 

Observations and applications of differential time delays in gravitational 
lenses are discussed in detail elsewhere in this book. Here some results which 
have implications for the millimetric part of the electromagnetic spectrum 
will be presented. More specifically, we will discuss a derivation of the dif- 
ferential time delay in the gravitational lens PKS1830-211 obtained from the 
saturated molecular absorption line of IIC0+(2-l). 

The gravitational lens system PKS1830-211 has been described in some 



detail above (Sect. 3.3). The background quasar is variable at radio wave- 
lengths, with an amplitude which increases at shorter wavelengths. This is 
due to the fact that the core, where the variability occurs, is a flat-spectrum 
source while the jet, which has a more or less constant flux, has a steep 
radio spectrum. It is presently unknown if PKS 1830-2 11 is variable at in- 
frared/optical wavelengths, although this is likely to be the case. 

Time delay measurements of the PKS 1830-2 11 system has also been done 



using long wavelength radio continuum |147| [g4|. In one case a single dish 
telescope was used and the two main lens components were not resolved |147|. 
The analysis had to be based on a compound light curve and the derived time 
delay of 44 ± 9 days should therefore be regarded as tentative. In the other 
case, the ATCA interferometer was used, with an angular resolution that 
did not fully resolve the NE and SW components |Q. Instead a model fit- 
ting procedure was used in order to obtain two separate light curves over 
an 18 months period. The resulting differential time delay is 26^5 days. Al- 
though the analysis is model dependent, this result represents a considerable 
improvement in the At estimate. 
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Time delay measurements using molecular absorption lines. As 

discussed in Sect. 3.3, the lens in the PKS1830-211 system was first detected 



through molecular absorption lines at a redshift = 0.88582 [151]. More 



than 16 different molecular species in 29 different transitions have so far been 
detected at millimeter wavelengths [151| [154] [56 . Two additional molecular 
species in three different transitions have been observed at cm wavelengths 

The millimeter transitions include three different isotopic variants: H^'^CO^, 
HC"'^^0+ and H^^CN. The mere detection of these lines shows that the main 
isotopic transitions of these molecules must be highly saturated. Despite this 
the absorption lines do not reach zero intensity (Fig. |l7|b). This can only be 
reconciled with an optical thick obscuration that do not completely cover the 
background continuum emission. In fact, from the ratio of the total continuum 
and the depth of saturated molecular absorption lines (such as HC0+(2-l), 
HCN(2-1), etc.), it was concluded that only the SW lens component is ob- 
scured by molecular gas and t hat the covering factor of this particular image 
is unity or close to unity ] [l5l| . A secondary weaker molecular absorption has 
now been found towards the NE component as well, separated in velocity by 
-147 kms-i [ p^ . 

Imaging of the HC0+(2-l) absorption line with the IRAM millimeter 
wave interferometer did not directly resolve the NE and SW components 
] |154| . This is due to the low declination of the source relative to the latitude 
of the interferometer, creating a synthesized beam elongated in approximately 
the same direction as the image separation. The continuum, however, is strong 
enough to allow self-calibration, making it possible to accurately track the 
phase center. The best angular resolution is achieved in right ascension (^ 
O'.'l) with a factor ^2 worse resolution in declination due to an elongated 
synthesized beam. At frequencies outside the absorption line, the phase center 
should fall on a line in between the NE and SW components. Assuming 
that the flux ratio NE/SW is similar to that derived for longer wavelengths 
(^ 1.3—1.4), the phase center should move towards positive RA at frequencies 
where the absorption occurs. If the covering factor is unity, Aa should be 
~ -1-0 '.'25. This is exactly the amount of shift observed for the saturated 
HC0+(2-l) line (Fig. [l7|a). A similar shift in the decl inat ion of the phase 
center can also be seen and concurs with these results ]154]. This result has 
also been confirmed through BI MA o bservations where the two continuum 
components have been separated ]140[. 

Due to this fortunate configuration of obscuring molecular gas in the lens- 
ing galaxy, the flux contributions from the NE and SW cores can easily be 
estimated using molecular absorption lines and a single dish telescope with 
low angular resolution. The 15m SEST telescope, which is used for the time 
delay monitoring presented here, has a HPBW of ^^50" at the observed fre- 
quency of the IIC0+(2-l) transition, much larger than the image separation 
of 0'.'97. Since molecular gas covers only the SW component, as shown by 
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the interferometric data, and the Une opacity is 3> 1 as seen from the rare 
isotopic Unes, the depth of the absorption line corresponds to the flux from 
the SW component only. The total continuum away from the absorption line 
corresponds to the sum of fluxes from the SW and NE components (Fig. |l8|) . 
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Fig. 17. Bottom: Spectrum of HCO+(2-l) at z = 0.886 towards PKS1830-211 
obtained with the TRAM interferometer. The main absorption line is seen around 
zero velocity. A secondary, weaker absorption line of HCO"'' (2-1) is seen at a velocity 
of —147 km/s relative to the main line. Top : The right ascension shift of the 
phase center of the continuum emission as a function of velocity. A negative shift 
means that the phase center moves towards the NE component, while a positive 
shift indicates a shift towards the SW component. Comparison with the absorption 
spectra shows that the main absorption component covers th e SW source, while 
the weaker secondary absorption covers the NE source. (From |154|. 



Monitoring of HCO+(2-l). Monitoring of the HC0+(2-l) absorption 
and the total continuum flux has been going on at the 15m SEST telescope 
since April 1996. Data can only be obtained between February and November 
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Fig. 18. Illustration of how the individual light curves of the NE and SW compo- 
nents can be derived fr om t he single dish millimetric observations. A 5 GHz radio 
image of PKS1830-211 [i3£] shows the two cores and the extended jet emission. At 
millimeter wavelengths only the cores contributes to the continuum emission. The 
right hand panel shows the saturated HCO^(2-l) spectrum. Since the molecular 
gas only covers the SW component and the line opacity is 3> 1, the depth of the 
absorption line corresponds to the flux from the SW component only. The total 
continuum away from the absorption line corresponds to the sum of the fluxes from 
the SW and NE components. 



due to Sun constraints (PKS1830-211 comes within 3° from the Sun). The 
light curves are shown in Fig. |l9[ Since only the total continuum (at the top in 
the figure) and the depth of the absorption line (at the bottom in the figure) 
are measured, the flux from the NE component (middle) is derived as the 
difference Tne = ?tot— Tabs, and therefore has a somewhat higher uncertainty. 
During the 1996-2000 campaigns a total of 144 usable observations have been 
obtained. The background quasar had a large outburst during 1998. The 
outburst shows a single peak in the total continuum, putting an upper limit to 
the differential time delay between the two cores (no double peak structure) . 
Separating the light curves for the two cores, however, one can clearly see a 
delayed response of the SW image relative to the NE (Fig. |l^). This is seen 
even more clearly when plotting the ratios of the NE and SW fluxes (Fig. po|) . 
This ratio shows the relative magnification of the two cores and should be 
constant in the absence of a time delay. The decrease in the fiux ratio during 
the 1998 outburst is a clear indication of the time delay between the NE and 
SW cores. 
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Fig. 19. Results from 5 years of monitoring of the HC0^(2-1) absorption at Zd = 
0.886 towards PKS1830-211. The top curve shows the measured total continuum 
flux away from the absorption line. Notice the large outburst during 1998. The 
bottom curve shows the depth of the HCO'''(2-l) line (shifted by —30 mK). This 
corresponds to the flux from the SW component. The middle curve shows the flux 
derived for the NE component. 
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Fig. 20. Same as Fig. |l^ but showing the flux ratio of the NE and SW components. 
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Monitoring results. The light curves shown in Figs. and have been 
analyzed using several different techniques. The problem is straightforward: 
correlate two unevenly sampled time series. However, the analysis is compli- 
cated by the fact that since the two time series are really copies of each other 
but shifted in time they are effectively sampled at different epochs. 

There are two main types of methods used for obtaining the time delay. 
One is to interpolate unobserved data points and then apply standard tech- 
niques for cross correlation. The other is to use the unevenly sampled time 
series and try to correlate neighboring points as good as possible. The former 
method is superior if the sampling rate is high, but this is usually not the 
case for astronomical data. The latter method has less precision but is the 
least unbiased way of obtaining the time delay. Both types of methods have 
been used for the molecular absorption line data on PKS1830-211. 



Data analysis. The observables at each epoch U are the total continuum 
flux So(ti) and the depth of the absorption line, corresponding to the flux 
from the SW component Ssw{ti)- The total continuum is the sum of the two 
image components, 

SoiU) ^ SiiU) + S2iU) , 

where subscript 1 refers to the NE component and 2 to the SW component. 
We know that the Si and ^2 fluxes are related as 

SiiU)^ ^^S2{t^ + At) , 

where /x is the magnification ratio and At the differential time delay between 
the two cores. Hence, 

SoiU)^ ^iS2{t^+At) + S2{U) . 

By shifting the observed ^2 values, multiplying it with a magnification ratio 
and adding the observed unshifted and unmagnified values we should recover 
the observed total flux at time t^. Since the observations consist of an unevenly 
sampled time series, with significant amount of noise, finding the true At and 
/X is a non-trivial exercise. 

The analysis of the light curves has been done using three methods. Two 
of them involves interpolating between the observed data points and con- 
struction of an evenly sampled time series. Due to the rather long inter- 
ruptions due to the Sun avoidance, interpolation only extends over periods 
between February and November. An elaborate interpolation scheme of un- 
evenly sampled data has been developed with the specific goal of resolving 



the time delay controversy of 0957-1-561 (Press et al. |117|). In the case of the 
molecular absorption data, however, a smoothing function with an effective 
resolution similar to the average data point separation at each epoch was 
applied. Using the smoothed time series, data points in between observed 
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epochs were linearly interpolated. The smoothing dampens the worst fluc- 
tuations while retaining the small scale structure in the time series. It also 
allows an easy assessment of the relative weights of observed and interpo- 
lated data. A complication, however, is that the data points are no longer 
completely independent. This is of some concern when deriving the reduced 
^ values. 



minimization : Minimization was done using both the time delay At 
and the magnification ratio /i as well as keeping the magnification ratio fixed 
or time dependent. The latter is due to a surprising realization that the 
magnification ratio might be variable, albeit on a much longer time scale 
than the time delay (cf. Fig. Using a parameterized means that the 
solution becomes cumbersome and slow. Instead we smoothed the flux ratio 
and fitted a third order polynomial. This parameterization of [liti) was used 
when solving for At. The result, together with results from the other analysis 
methods, is shown in Fig. Ell and gives & At — 21 days. 



Cross correlation : Edelson & Krolik |46| developed a discrete cross cor- 
relation method specifically aimed for reverberation mapping of AGNs that 
can be used for time delays in gravitational lensing. The method optimizes 
the binning of data points rather than the interpolation, as in the method of 
Press et al. |117|. The method requires a fairly well sampled data set to start 
with in order to retain a sufficiently good temporal resolution. The sampling 
rate for molecular absorption line data in PKS1830-211 is not dense enough 
to use the Edelson & Krolik method. Instead cross correlation was done on 
the same smoothed and interpolated data set as the minimization. The 
cross correlation coefficient is defined as Tab = •Sab/('^a Sb), where the co- 
variances Sa and Sb are defined in the usual manner (cf. [|l5j). As with the 
minimization, the data points are not entirely independent due to the 
smoothing and interpolation and the variances are only approximately true. 
The result gives At = 25 days, with a rather broad maximum for the cross 
correlation coefficient (Fig. |2l| ). 



Minimum dispersion (the Pelt method) : A simple and robust tech- 
niq ue for ana lyzing unevenly sampled time series was presented by Pelt et 
al. |110| llll]. They successfully applied it to the lens system 0957-1-561. The 
strength of the method is that interpolation or smoothing are not needed, 
leaving the errors for each data point independent. The method is a form 
of cross correlation where a given data point is correlated with a data point 
which is temporally its closest neighbor. The method is illustrated in Fig. ^2], 
where the round and square markers in the two top rows represent the two 
photometric data sets obtained from a two-component gravitational lens. 
When correlating the time series, one of them is shifted in time, as the square 
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markers in the middle rows. Projecting both the unshifted (round) and the 
shifted (square) time series to a common array (bottom row), correlation is 
done between those data points which are from different time series and clos- 
est to each other. These points are connected by arcs in the figure. It is easy 
to include the effects of different magnifications for the lensed components as 
well as time delays in systems with two or more lenses [ |112| | . 

The results are undeniably noisier than for the interpolated data sets. This 
can be seen in Fig. where the Pelt dispersion method has been applied 
to both raw and interpolated data. The best fit is for a time delay At = 28 
days, with the NE component leading. 



Error analysis : The errors associated with the light curves are a combina- 
tion of noise in the data points and systematic errors. The latter can originate 
in the instrument, in the modeling necessary for separation of the lensed com- 
ponents (as in the case of long wavelength radio observations), assumptions 
made about the lensing system, etc. When interpreting the time delay in 
terms of a Hubble constant, the largest systematic error comes from model- 
ing of the gravitational potential (see Chapter X). Noise in the data comes 
from imprecise measurements but can also originate in secondary variability 
such as microlensing and interstellar scintillation. The latter is applicable at 
long radio wavelengths. Microlensing may be of importance even for gravita- 
tional lenses observed at radio wavelengths (cf. [^). 

In order to assess the significance of correlations found in the light curve of 
gravitationally lensed images it is customary to derive the confidence limits 
through Monte Carlo simulations and bootstrap techniques (cf. [|l)). The 
results often shows non-Gaussian distributions and confidence levels are set 
by finding the range of delays and magnification ratios inside which a given 
amount (say 95%) of the simulations lie. This gives a better estimate of the 
true confidence level than simply fitting a Gaussian to the distribution. Doing 
this for the molecular absorption line data in PKS 1830-2 11 gives a time delay 
of At = 28^ days, with the NE component leading. 

In Fig. g4|the light curve of the SW component in PKS1830-211 has been 
shifted by —28 days and multiplied by a magnification ratio ^. In the upper 
panel a constant ratio of ^ = 1.3 was used, while in the lower panel a time 
dependent magnification ratio was used. The use of different parameteriza- 
tions of the magnification ratios do not change the derived time delay, but the 
time dependent form provides a better fit of the two light curves. The reason 
for the slow change in magnification ratio is presently unclear. It may have 
implications for the use of molecular absorption lines as a probe of the time 
delay, but since the time scale for the change of the magnification appears 
to be much longer than the time delay, it is likely to be of small importance 
when correlating light curves for each period (i.e. 9 months). 
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Fig. 21. The time delay for PKS1830-211 derived from molecular absorption lines 

and using cross- correlation, minimization and the Pelt Minimum Dispersion 
method. All three methods in this example use the interpolated data set. 
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Fig. 22. Illustration of the Pelt Minimum Dispersion method used to determine 
time delays in gravitational lenses. This method was one of the methods used 
for deriving the time delay in PKS1830-211 from molecular absorption lines. This 

particular case shows a two-image lens (A and B), where the respective light curves 
are sampled at irregular intervals. In the middle section the B light curve is shifted 
by At. By projecting the resulting data points to a common array (C), nearest 

neighbors of different light curves are correlated (arcs). A weight, depending on the 
time difference between the points used in the correlation, can be applied. 
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Fig. 23. Results for the time delay of molecular absorption lines in the PKS1830- 
211 gravitational lens system using the Pelt Minimum Dispersion method. Both 
raw and interpolated data sets are shown (left and right, respectively). Also shown 
are the results for different treatments of the magnification ratio. In the top two 
rows, the magnification ratio is included in the fit, while in the two bottom rows 
the magnification ratio is predetermined by fitting either a second or third degree 
polynomial to the observed magnification ratio. 



6 LENS MODELS FOR PKS1830-211 

In order to use the differential time delay to derive a value for the Hubble 
constant, a lens model has to be fitted to the observed data. This is not 
a trivial exercise in most cases and this is particularly true for PKS1830- 
211. Due to its location at Galactic longitude / = 12.2° and latitude b = 
—5.7°, PKS1830-211 suffers considerable Galactic extinction. In addition, 
the molecular gas seen in absorption towards the SW component contributes 
significant obscuration for at least this image. Early attempts to identify the 
radio source PKS1830-211 with an optical counterpart were all unsuccessful 
[ [139| pOf . It was only with the advent of sensitive infrared imaging and 
spectroscopical capabilities that progress could be made. The NE image was 
positively identified using K-band imaging at Keck and the ESO NTT pof . 
While the redshift of the lens, Zd = 0.886, had been derived using molecular 



absorption lines |151], the redshift of the source was obtained from near- 
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Fig. 24. Tlie light curve of the SW component in PKS1830-211 (black) shifted by 
At = —28 days and multiplied by a magnification ratio ^. The light curve for the 
NE component is shown in grey. In the upper panel a constant magnification ratio 
of /J = 1.3 was used, while in the lower panel a time dependent magnification ratio 
was used (cf. Fig. Ei). 



infrared spectroscopy ||9j. The redshift was found to be Zs = 2.507. Imaging 
with the HST WFPC2 and NICMOS allowed identification of both the NE 
and SW image [0. In addition, an object which might be the lensing galaxy 
was detected (designated as G). Its exact center position remains uncertain 
due to the presence of a point source ~ 190 mas away. The nature of the point 
source remains unknown but could possibly be a Galactic star and thus of 
no importance for the lens model. In a recent paper by Courbin et al. 
combined images from the HST and Gemini-North telescopes show what 
might actually be the lensing galaxy at Zd = 0.889. The lens has two spiral 
arms, as expected from the molecular absorption data. One spiral arm crosses 
the SW image of the QSO. The center of the spiral is, however, significantly 
offset from the line joining the NE and SW images. Based on symmetry 
arguments, the center of the lensing galaxy is believed to be in the proximity 
of this line joining the images (Fig p5| ). 

All the necessary ingredients for a detailed lens model are thus in place, ex- 
cept for two remaining uncertainties: the exact position of the lensing galaxy 
and the possible double-lens nature of the system (cf. Sect. 3.3). 
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Fig. 25. Illustration of the positions of various components in the PKS1830-211 
gravitational lens system. The two cores are marked by QSO (NE) and QSO (SW). 
The putative location of the lens galaxy, derived by Lehar et al. IstI, is marked by 
'G' (filled circle) and a point source of unknown origin is marked by 'P' (star). A new 
possible lens location (from ^]) is marked by 'S' (filled triangle). The location of 
the center of the lens derived from the lens model here, using the Lucy rectification 



scheme, is marked by 'Lucy' (open circle), the lens center derived by Nair et al. [102 
by 'NRR93' (open diamond) and the lens center derived by Kochanek & Narayan 
]7E| ] by 'KN92' (open rectangle). The dashed line joining the two QSO images is 
only to guide the eye. 



6.1 Early models 

There have been several attempts to model the lens system PKS1830-211. 
The system was first detected at radio wavelengths, where it is a prominent 
southern radio source. The morphology of the system was found to be that 
of a double, while the radio spectrum is typical for a compact flat-spectrum 



source. This led Rao & Subrahmanyan |120 to first suggest that PKS1830- 
211 is a gravitationally lensed system. Based only on the radio images and 
their polarization properties, obtained with the Very Large Array (VLA) at 



5 and 15 GHz, Subrahmanyan et al. |139] constructed a lens model which 



is not much different from later ones based on more detailed data. In order 
to reconstruct the extended radio structure, Subrahmanyan et al. modeled 
the source as an one-sided core-jet structure. To get lensed images with a 
morphology similar to the observed one, they also had to include a 'knot' 
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in the jet. Based on their lens model Subrahmanyan et al. predicted a time 
delay of 27/ij~QQ days (using the now known redshift of the source and the 
lens). 



Nair et al. |102| modeled the PKS1830-211 system using improved radio 
interferometry data (cf. [^). The method was similar to that of Subrah- 
manyan et al. |139| in that the source structure was built up in a piece- 
meal manner in order to fit various observed features. With this type of 
method one can emphasize the influence of small and weak features which 
may carry a small weight in an inversion scheme based on minimization of 
model— observed results, but may nevertheless carry important information 
on the lensing scenario. In the PKS1830-211 system such a weak radio feature, 
labeled E, was used by Nair et al. to constrain the lens model. This feature 
might be a third demagnified image of the core. However, since the flux of 
the E component is less than one percent of the peak value, its signiflcance 
in terms of flux is small unless the dynamic range of the interferometry maps 
is very good. Nair et al. found that an elliptical potential with the radio core 
located close to the inner edge of the radial caustic gave a good flt to the ob- 



served morphology. As in the previous model by Subrahmanyan et al. |139|, 
it was necessary to include a 'knot' feature in the source distribution. The 
jet needed to be bent and cross the tangential caustic. The model gives a 
good fit to the observed system and places the lens galaxy close to, but not 
coinciding with, the possible lens position observed by Lehar et al. [|^. The 
estimated time delay, using the known source and lens redshifts, is 
days. 

A difficulty with extended lensed images is that any inversion must solve 
simultaneously for both the lens configuration and the source structure. This 
type of inversion problem can be seen as 

lobsiC)^ I HC)K{C-C)dC , (7) 

where both the source distribution V'(C) the kernel K (here representing 
the lensing potential) are unknown. This type of problem is generally unsolv- 
able. In the case of lensing, however, one can use the knowledge that when 
the lensed image contains multiple distorted components of the background 
object these must arise from a common source. Furthermore, it is known 
that surface brightness is conserved. These 'priors' constrain the problem and 
permit the simultaneous solution of both the structure of the source and the 
properties of the lensing potential. This type of inversion problems for gravi- 
tational lenses has been developed extensively by Kochanek, Wallington and 
collaborators in several papers (cf. |^ jT^ |148|). In particular, Kochanek 



& Narayan ||79|| developed an inversion method based on the CLEAN rou- 
tine (cf. psf) used in radio interferometry data reduction and applied it to 
PKS1830-211. This method takes into consideration the effects of the finite 
resolution when attempting to invert the lens model and is thereby able to 
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better distinguish the best lens model. The LensClean method of Kochanek 
& Narayan has produced the hitherto most reliable model for the PKS1830- 
211 system, but due to the finite resolution, the inversion was done on radio 
data with rather low angular resolution but with good signal-to-noise, it is 
not likely to represent the final model. 

Lehar et al. |Q modeled the PKS1830-211 system using a singular isother- 
mal elliptical mass distribution as well as with two singular iso thermal spheres 
representing the lens galaxy and the source G2 (cf. Sect. 3^). In both these 
cases they fixed the lens at the position of G, with a positional uncertainty 
of 80 mas. The extended radio emission was not used to constrain the lens 
model. Lehar et al. noted the strong dependence of the location of the lens 
galaxy and the Hubble constant derived from differential time delay measure- 
ments. 
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Fig. 26. A 15 GHz radio image of PKS1830-211 (from ||l39|). The center of the 
coordinate system is arbitrary. All positions in the text and in Table ^ are relative 
the NE image. The critical lines of the best fit lens model, shown in Fig. ^ are 
shown for comparison. 
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Fig. 27. The lensed images obtained from inversion of the lens equation using Lucy 
rectification and Simulated Annealing (as described in the text). In this particular 
solution, the lens is fixed at the position of the observed (putative) lens center G 
(see Fig. Ga). The critical lines from the lens model are marked. 



6.2 A new lens model of PKS1830-211 

The situation is rather unsatisfactory concerning the various solutions to 
the le nsin g configuration characterizing PKS 1830-2 11. Three different models 
(]79|] 1 102 1 |8^]) give three different positions for the lens, with corresponding 
differences in the value of the Hubble constant for a given time delay. The 
situation is summarized in Fig. ^ The location of the cores in the long 
wavelength radio data used by Kochanek & Narayan differ from that used 
by others (using data from shorter wavelengths). Their lens center relative 
to the other features shown in Fig. ^ is therefore somewhat uncertain. 

The fact that surface brightness is conserved can be used in modeling a 
lensing configuration if both the observations and the code have infinitely 
good resolution. In the opposite extreme, if the source remains unresolved, 
one can calculate the magnification for a point source and apply a smoothing 
function (convolution) representing the observational transfer function (e.g. 
atmosphere, telescope, imaging array). A more difficult situation arises when 
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Fig. 28. The source distribution derived from the best fit lens model, which gives 
the image seen in Fig. The caustic structure is sho-wn. 



an extended source distribution is partially resolved by the observer. A given 
resolution element will represent different areas of the source in a rather 
complicated manner and the conservation of surface brightness ceases to be a 
good prior for an observed lens system. In their LensClean method Kochanek 
& Narayan Jzof solves this situation by representing the source emissivity 
distribution by ^-functions, mapped through the lens configuration with the 
corresponding magnification 'turned on' and then smoothed by a restoring 
beam. 

An alternative way to solve the inversion problem as stated above is pre- 
sented here. A more thorough description of the method will be published in 
Wiklind (2002). The concept is similar to the LensClean method of Kochanek 
& Narayan [Q. However, instead of introducing CLEAN components, as S- 
functions to represent the source distribution, this method starts with a none- 
zero smooth source distribution and applies the Lucy rectification method [p5| 
to constrain the source emissivity, given the observed lensed images and for 
a given lens model. The Lucy rectification scheme has been used extensively 
to deconvolve images obtained with the Hubble Space Telescope before the 
corrective optics was installed. It has also been used to deconvolve the molec- 
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ular gas distribution in galaxies observed with single dish telescopes of rather 



poor angular resolution |156| [158|. 



The method has an outer loop, which controls the lensing parameters, 
and an inner loop which solves for the best source distribution given the lens 
configuration. In the inner loop the source distribution is mapped through 
the lens and the resulting image is compared with the observed one. The 
source distribution is adjusted according to the Lucy method (see below) 
and the process is repeated. When the inner loop has converged, the source 
distribution is mapped through the lens a final time and the resulting image 
gives a goodness-of-fit. The lensing configuration is then modified in the outer 
loop and the process is repeated. The Lucy method is used only in the inner 
loop, each time starting with a perfectly smooth source distribution. The 
observed images are represented by 

Iobs{xi,X2) ^ J J tp{yi,y2)K{xi,x2\yi,y2)dyidy2 , (8) 

where ip{yi, j/2) is the true source distribution and K{xi,X2\yi,y2) is the ker- 
nel representing both the gravitational lens and the finite angular resolution 
of the observation. The kernel is here written as a conditional probability 
function: the likelihood of (xi,X2) given (2/1,2/2)- With this formulation we 
can use the Lucy method in a straightforward manner. The idea being that 
an approximation to the true source distribution is 

r+\yi,y2)^riyi,y2) J j ^-f^^^ K{x^,x2\yuy2)dx,dx2. (9) 

The conditional probability function K contains the likelihood of an image 
emissivity at position (a;i,a;2) given a source emissivity at (2/1,2/2) and the 
action of a restoring beam (i.e. a finite angular resolution). The simplest 
(although not entirely correct) restoring beam is a Gaussian with a HPBW 
similar to the angular resolution of the observations. 

Even starting with a constant source emissivity distribution V'", the Lucy 
rectification converges very rapidly to a specific source distribution. Unfortu- 
nately, there is no good criteria for determining when to stop the rectification 
(cf. [|95[ ). In this particular application the inner loop was stopped after 5 it- 
erations of the Lucy rectification. The outer loop consists only of changing the 
lensing parameters. Several different methods can be employed for this, the 
most efficient for this application being the Simplex method. However, there 
is a risk that this method gets stuck in a local minimum and great care has to 
be taken to ensure that a global minimum has really been reached. This in- 
volves repeatedly restarting the Simplex method with parameters offset from 
the ones giving a (local?) minimum in the xt. An alternative method that 
circumvents this, but that is computationally more expensive, is Simulated 



Annealing (cf. |99[| |15£|). This latter method was used in this application. 
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The lens was modeled as a non-singular elliptical mass distribution 



k{xi,X2) = KO {^^i + ^ + s"^ j , (10) 

where q is the projected axis ratio, s is the core radius and the surface density 
profile is set to 7 = 1/2, representing an isothermal mass distribution. The 
deflection angle and magnification was calculated using the code developed 
by Barkana This code can handle surface density profiles with 7 ^ 1/2, 
but here the modeling is restricted to the isothermal case. The density profile 
is, however, a very important parameter when deriving the Hubble constant 



using differential time delays (cf. |159| ||8^). In all there are six lens param- 



eters that were fitted: the center position of the lens, the position angle, the 
velocity dispersion, the ellipticity, and the core radius. 

Applyin g thi s method to the 15 GHz radio image of PKS1830-211 shown 



in Fig. |2^ 139 , a best fit lens model is achieved with the parameters as 
tabulated in Table ||. Also listed in Table || are the results when keeping the 
lens position fixed at the coordinates of the putative infrared lens center G. 
The resulting image distribution for case (b) is seen in Fig. ^ together with 
the critical lines of the lens. The corresponding source distribution is seen in 
Fig. |2^ together with the caustic structure. 

The two solutions presented in Table |^ are very similar to each other, yet 
they give quite different values to the Hubble constant. Using the lens model 
of Nair et al. |102|, the Hubble constant becomes Hq = 59^^^ km/s/Mpc. 



These differences are mainly due to the different locations of the lens center 
and introduces a large uncertainty in the correct value of Hq. 

The lens model presented here will be further refined and the results 
should be regarded as tentative. However, unless the position of the lens can 
be derived more accurately, the value of the Hubble constant will remain 
uncertain. As mentioned above, the shape of the density profile is also a 
source of uncertainty for a more exact derivation of Hq. This uncertainty is 
largest for exponents 7 < 1/2 p3|. 



7 FUTURE PROSPECTS 

Existing millimeter and submillimeter telescopes use gravitational lenses more 
as an aid to the study of distant objects, rather than being an aid to the 
study of gravitational lensing as such. Nevertheless, some information about 
the content of the interstellar medium in both lenses and sources have been 
obtained, and the molecular absorption lines seen towards PKS 1830-2 11 have 
been used to measure the differential time delay in this particular system. 

This situation will change dramatically when planned telescopes at mil- 
limeter and submillimeter wavelengths become available. Increased sensitivity 
and angular resolution will make this wavelength regime very important for 
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Table 3. Lens parameters 





6 free parameters 


4 free parameters 
(lens fixed) 


lens center" 


-0.5008 -0.5205 


-0.5010 -0.4450'' 
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(a) Relative to the NE component. 

(6) Fixed at possible lens center G |8^ ]. 

(c) With At = 281:1 days. 



studies of gravitational lenses as a phenomenon of their own. The most ob- 
vious advantage is that obscuration effects will be completely absent. The 
effects of microlensing will also be absent or at least minimal. The use of 
flux ratios of lensed components for constraining parameters when modeling 
lensed have fallen out of favor due to differential extinction and microlens- 
ing effects, but will be usable when the new submm/mm instruments are 
available. 

Existing submillimeter and millimeter facilities include single dish tele- 
scopes, such as the IRAM 30m telescope on Pico Veleta in Spain, the JCMT 
15m telescope on Mauna Kea, and the SEST 15m telescope on La Silla in 
Chile. Two 10m size dishes, aimed primarily for submillimeter wavelengths, 
include the CSO on Mauna Kea and the HHT on Mount Graham in Arizona. 
These telescopes use both heterodyne receivers for spectral line observations 
and bolometer type array cameras for continuum observations. The sensi- 
tivity depends largely on the quality of the site and the instrumentation. 
The angular resolution, however, is determined by the diffraction limit of 
the telescopes. At A = 1mm, the angular resolution is limited to 10" — 25". 
This constitutes the largest l imit ation to the study of gravitational lenses. 
For number counts (see Sect. 4.2) the lack of angular resolution means that 
with only slightly more sensitive receivers, confusion will become a major 
limitation (cL |l65| ) 

A few interferometers operating at millimeter wavelengths exist. The 
IRAM Plateau de Bure interferometer in France consists of five (soon to 
be six) 15m telescopes, and represents the largest collecting area today. The 
OVRO interferometer consists of six 10m telescopes, while BIMA consists 
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of eight 6ni telescopes. In Japan the Nobeyama interferometer consists of 

six lOm telescopes. All of these facilities operate at A = 3 — 1mm. The an- 
gular resolution reached is typically around 1" or slightly better. However, 
sensitivity becomes a serious limitation at the longest baselines and highest 
angular resolutions. Also, the Australian Telescope Compact Array (ATCA) 
has recently been upgraded to work at 3mm with five of its 22m elements. 



10' 




Fig. 29. A comparison of the wavelength coverage and the angular resolution of 
existing and planned instruments. The VLA refers to the Extended VLA. With 
ALMA, the same angular resolution will be reached from cm wavelengths into the 
optical. 



7.1 Future instruments 

Single dish telescopes. A few single dish submm/ mm telescopes are under 
construction, or in advanced planning. These include the recently commis- 
sioned Green Bank Telescope (GBT), which has a 90m unobstructed dish and 
will reach A3mm when fully operational. The Large Millimeter Telescope is 
a 50m telescope built as a collaboration between INAOE and the Univer- 
sity of Massachusetts in Amherst. This telescope will be placed on Sierra 
Negra in Mexico and operate at a wavelength of 1-4 millimeter. Two single 
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dish telescopes to be situated close to the ALMA site (see below) are under 
construction: APEX is a 12m, single dish telescope which will be placed on 
Chajnantor at an altitude of 5000m. This telescope will operate into the THz 
regime, i.e. A ^ 300^m. ASTE is a Japanese 10m dish, to be placed at Pampa 
La Bola, a few kilometers away from the ALMA site. The ASTE will also 
operate at submm to THz frequencies. These new single dish telescopes will 
explore distant objects, including gravitational lenses, with somewhat better 
angular resolution than existing telescopes. Nevertheless, even at the highest 
frequencies it will only reach an angular resolution of 6". This is insufficient 
for detailed studies of gravitationally lensed systems. 



The Atacama Large Millimeter Array. A major step in submm/mm 
wave instruments will be the joint European-US project, with Japanese in- 
volvement as well, of building a large millimeter and submillimeter interfer- 
ometer at an altitude of 5000m on Chajnantor in Chile. This instrument, 
with the acronym ALMA, will consist of 64 x 12m telescopes, each with a 
surface accuracy of at least 20/xm. The total collecting area will be 7238 m^ , 
which is an order of magnitude greater than the largest existing instrument 
today. The longest baseline will be 10 — 12 km, leading to an angular resolu- 
tion surpassing that of the Hubble Space Telescope. A rough estimate of the 
angular resolution power is 0'.'2 Amm/Lkm, where A is the wavelength in mil- 
limeters and L is the baseline length in kilometers. ALMA will easily reach 
O'.'l resolution and, at the highest frequencies and longest baselines, O'.'Ol. A 
comparison of the wavelength coverage and the projected angular resolution 
for existing and planned telescopes is shown in Fig. |2^. 

At the same time ALMA will increase the sensitivity over existing instru- 
ments by at least two orders of magnitude. The noise rms level reached with 
an interferometer consisting of N array elements, each of effective area Acs 
and with an integration time of tint seconds over an effective bandwidth of 
i? Hz, is expressed as 

AS=^^=M^^= , (11) 
^N{N-l)BU^t 

where the term e'^ represents the damping by the atmosphere at frequency 
1^. The system temperature, T^yg, is the total noise received by the telescope 
(including ground pick-up). From this expression it is clear that large indi- 
vidual antenna sizes, a large number of array elements, a large bandwidth, a 
low system temperature and as little damping from the atmosphere are es- 
sential for a sensitive interferometer. These criteria can be fulfilled by using 
state-of-the-art receiver systems and putting them at a high and dry site such 
as Chajnantor in the Andes. 

The excellent atmospheric conditions at Chajnantor is estimated to give a 
total system temperature (including the atmospheric damping) of ^^230 K at 
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a wavelength of 0.8mm. Using a 16 GHz wide backend and one hour of inte- 
gration, the 5(7 detection hmit will be 100 /xJy. At a wavelength of 460 /im the 
corresponding 5(t detection limit after one hour is 1 Jy. Although considerably 
worse, this wavelength range is usually not reachable at all with ground-based 
facilities. These values are representative for unresolved sources. If the source 
is resolved at the longest baselines, the sensitivity decreases. 




Fig. 30. A simulated 9 arcmin^ deep field observed with ALMA at a wavelength of 
850/im. All sources above the 5a noise rms hmit of 0.15 mjy are shown. The source 
population was modeled to fit the observed cumulative source counts, extrapolated 
to 0.1 mJy using a more shallow powerlaw than what can be fitted in the range 
1-10 mJy. Left: Unlensed field. Approximately 210 sources are 'detected'. Right: 
The same field as in the left image but with an intervening rich galaxy cluster at 
redshift z = 0.3. 



7.2 Weak lensing at submillimeter wavelengths 

With its superb sensitivity and angular resolution, ALMA can successfully 
be used to study weak lensing by intermediate redshift galaxy clusters. As- 
suming a constant co-moving volume density of dusty high redshift objects, 
the constant sensitivity to dust continuum emission from z « 1 to z « 10 
(Sect. 4.2) favors the detection of the highest redshift objects. This selection 
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bias is unique for the submm/mm wavelength regime. Furthermore, a reason- 
ably well-sampled ww-plane, resulting from a large number of interferometer 
elements, will have a point-spread function (PSF) which is well behaved in 
comparison to existing wide-field CCDs. And, finally, the dust distribution in 
galaxies tend to be symmetric and centrally concentrated, leading to a simple 
geometry of the lensed sources. The latter issue is, of course, not well con- 
strained by existing data. Although dust continuum emission will be easier to 
detect at large distances than line emission (cf. ) , the high sensitivity and 
broad instantaneous bandwidth of ALMA will allow CO lines to be observed 
for many of the lensed sources. This enables a direct way of deriving the red- 
shift distribution of the source population. Furthermore, observations of CO 
lines in the intervening galaxies will allow a determination of the dynamical 
mass of the lenses. 

ALMA will, however, have a relatively limited instantaneous field of view. 
This is particularly aggravating at high frequencies. The field of view is lim- 
ited by the size of the individual telescopes and will be similar in size to 
that of existing single dish telescopes. This means that the time consuming 
process of mosaicing is necessary. The key parameters for the effectiveness or 
speed of this new instrument are the sensitivity AS, the angular resolution 
and the primary beam The time required to survey a given area 

A to a flux density limit S is then: t (AS/ S)^{A/Afo^) (cf. |^ 0). 

An example of what can be achieved with ALMA is shown in Fig. 
The left panel shows a simulated 9 arcmin^ field observed with ALMA to a 
5(7 noise rms limit of 0.15 mJy. The wavelength is 0.8mm and was chosen to 
maximize the number of detected dust continuum sources while minimizing 
the observing time (by making the field-of-view as large as possible). The 
field of view per pointing is ~0.07 arcmin^. As shown in Sect. it takes 
somewhat less than one hour of telescope time to reach this 5cr detection limit. 
To cover 9 arcmin^ thus requires a total of 130 hours. In this example more 
than 200 sources will be detected (left panel) and their distortion by a rich 
cluster at z = 0.3 will be clearly detectable (right panel). In this example 
a population of dusty galaxies was assumed to have a constant co-moving 
volume density between z — 1 — 7 and undergo pure luminosity evolution. 
The observed flux densities were fitted to the number counts derived from 



SCUBA and MAMBO observations (Sect. The median luminosity of 

the 'detected' galaxies is 3 x 10^^ L0 and their redshift distribution is more 
or less flat between z = 2 to z = 7. The cumulative source counts was not, 
however, extrapolated to weaker flux densities using the same powerlaw as 
applicable between 10-1 mJy (cf. Fig. |8[ but with a more modest slope. This 
ensures that the number of sources seen in the right panel of Fig. ^ is not 
an overestimate. 

In the same manner as lensing by intermediate redshift clusters is presently 
being used to observe dust continuum sources at flux density levels otherwise 
not reachable with existing instruments, ALMA will reach flux densities ap- 
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preaching /uJy levels when observing lensed sources. For high redshift galax- 
ies, this flux density level corresponds to luminosities around 10^ L©, i.e. 
dwarf galaxies. 

8 SUMMARY 

The study of gas and dust at high redshift is important for several reasons. 
It gives us an unbiased view of star formation activity in obscured objects 
and it tells the story of the chemical evolution and star formation history in 
galaxies through the amount of processed gas (and dust) it contains. With 
today's millimeter and submillimeter facilities, this research area has used 
gravitational Icnsing mostly as a tool to boost the sensitivity. This is evident 
through the preponderance of gravitationally lensed objects among those 
which have been detected at 2; > 2 in the lines of the CO molecule. It is 
also evident in the use of Icnsing magnification by galaxy clusters in order to 
reach faint submm/mm continuum sources. There are, however, a few cases 
where millimeter lines have been directly involved in understanding lensing 
configurations. The best example of this is the highly obscured PKS1830- 
211, where the lens was identified through molecular absorption lines and 
where these lines give a velocity dispersion measure by originating in two 
different regions of the lens. The molecular absorption lines in this system 
have also been used to derive the differential time delay between the two main 
components, the main objective being to determine the Hubble constant, but 
also adding to the constraints in modeling this particular lens system. 

With future millimeter and submillimeter instruments, such as ALMA, 
coming on-line, the situation is likely to change drastically. The sensitivity of 
ALMA will be such that it does not need the extra magnification from lensing 
to observe very distant objects. Instead it will be used to study the lensing 
itself. The more or less constant sensitivity to dust emission over a redshift 
range stretching from z « 1 to z « 10 means that the likelihood for strong 
lensing of dust continuum detected sources is much larger than for optically 
selected sources. ALMA will therefore discover many more lenses and allow 
a direct assessment of cosmological parameters through lens statistics. Weak 
lensing will also be an area where ALMA can successfully contribute. Again, 
the high sensitivity to dust emission out to very high redshifts, combined with 
an angular resolution < O'.'l, and a more beneficial 'PSF' will make ALMA 
more efficient for probing the potential of galaxy clusters than present day op- 
tical/IR telescopes. In addition we will be able to study both the sources and 
the lenses themselves, free of obscuration and extinction corrections, derive 
rotation curves for the lenses, their orientation and, thus, greatly constrain 
lens models. 
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